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SUMMARY 


i 


A parallel  shock  wave  structure  is  modelled  as  a region  of  Inter- 
penetrating streams  of  cold  unshocked  and  hot  shocked  ions.  Our  linear 
dispersion  theory  predicts  that  unstable  whistlers  can  stand  and  therefore 
grow  to  large  amplitude  at  the  leading  edge  of  weak  and  Intermediate  strength 
shocks. 

Magnetosonic  instabilities  and  their  role  in  the  structure  of  perpen- 
dicular shock  waves  were  also  studied.  Since  these  instabilities  are  not 
stationary  in  the  shock  layer,  they  cannot  play  a principal  role  in  the 
structure. 

The  linear  dispersion  analysis  of  parallel  shock  structures  has  been 
used  to  estimate  the  fluxes  of  energetic  electrons  injected  into  the  mag- 
netosphere due  to  multiple  nuclear  bursts.  It  has  been  found  that  fluxes 
of  20  kev  electrons  two  orders  of  magnitude  larger  than  those  due  to  mag- 
netospheric  substorms  might  be  expected  to  occur  for  times  as  long  as  30  s. 
These  could  play  a significant  role  in  charging  synchronous  orbit  spacecraft. 
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A.  INTRODUCTION 

During  a high  altitude  nuclear  burst,  debris  particles  are  ejected 
from  the  fireball  by  essentially  two  processes.  The  first  process,  which 
is  collision  dominated,  involves  charge  exchange  interactions  between  the 
expanding  debris  ions  and  background  air  neutrals.  This  results  in  the 
ejection  of  debris  neutrals  from  the  fireball;  it  is  unlikely  that  this 
neutral  debris  ejection  can  be  directly  connected  with  any  mechanism  leading 
to  electron  injection  into  the  magnetosphere.  The  second  process,  collision- 
less in  nature  and  involving  turbulent  coupling  between  the  ejected  debris 
and  background  air  plasmas,  does,  however,  appear  to  bring  about  electron 
streaming  into  the  Van  Allen  belts  and  has  therefore  been  of  interest  to  us. 

For  example,  in  the  presence  of  the  compressed  magnetic  field  which 
can  penetrate  the  debris  bubble,  the  plasma  turbulence  takes  the  form  of 
large  amplitude  whistler  waves.  These  waves  can  suitably  couple  the  back- 
ground air  plasma  with  the  expanding  debris-air  plasma.  More  and  more  hot 
electrons  are  thus  acquired  and  some  have  sufficiently  large  velocities  to 
escape  along  the  distended  field  lines.  Moreover,  the  nonlinear  wave- 
resonant  particle  interactions  should  produce  anomalous  resistivity  whose 
scale  length  determines  the  extent  of  magnetic  field  penetration  into  the 
bubble  which,  in  turn,  determines  the  rate  of  escape  of  6 and  plasma  electrons. 

Superalfvgnic  debris  plasma  can  also  escape  directly  through  loss  cone 
exits  in  the  debris  bubble  since,  in  general,  the  cylindrical  axis  of  the 
bomb  casing  is  not  initially  aligned  with  the  direction  of  the  geomagnetic 
field.  The  more  perpendicular  the  cylindrical  axis  is  to  the  field  lines, 
the  greater  the  number  of  such  escape  particles.  Their  superalfvdnic 
velocities  suggest  the  formation  of  parallel  collisionless  shock  waves 
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^tshock.  II  £ = geomagnetic  field).  Our  studies  (Refs.  1,  2)  reveal  that 
such  shock  fronts  are  structured  by  turbulent  whistler  modes  which  couple 
the  incoming  background  air  plasma  to  the  shocked  debris  plasma.  Air 
plasma  can  therefore  be  picked  up  by  the  loss  cone  debris  and  deposited  into 
the  Van  Allen  belts. 

Evidently,  the  presence  of  a magnetic  field  has  a significant  effect 
on  shock  wave  structure.  First,  gradients  in  the  magnetic  field  give  rise 
to  electron  currents  that  can  drive  ion  acoustic  waves  unstable  and  increase 
the  effective  collision  frequency  (Refs.  3,  4)  (this  should  also  dictate 
the  penetration  depth  of  the  compressed  magnetic  field  into  the  debris 
bubble,  so  that  the  rate  of  escape  of  debris  and  air  electrons  could  be 
profoundly  affected).  Second,  when  propagation  is  perpendicular  to  the 
magnetic  field,  the  magnetic  field  can  inhibit  the  electrons  from  shorting 
out  ion  plasma  oscillations  for  wavelengths  long  compared  with  the  electron 
gyroradius  (Refs.  5 to  10);  for  wavelengths  short  compared  with  the  electron 
gyroradius  (Refs.  11,  12),  interactions  between  the  ion  beam  mode  and  the 
electron  Bernstein  modes  generate  instabilities  which  are,  however,  stabilized 
by  electron  heating,  resonance  broadening,  or  ion  trapping.  Third,  for 
oblique  or  parallel  propagation,  interactions  of  whistler  waves  with  ion 
acoustic  beam  modes  (Ref.  13)*  or  with  ion-cyclotron  beam  modes  (Refs.  1,  2) 
are  likely  to  be  important  and  the  existence  of  whistlers  assuredly  depends 
upon  the  presence  of  a magnetic  field. 

During  the  two  year  period  of  this  contract,  we  have  investigated 

several  streaming  plasma  instabilities  and  their  roles  in  the  structures  of 

*The  s called  "modified"  two  stream  instability  described  in  Refs.  7 to  9 
is  operative  for  propagation  angles  slightly  less  than  90°,  so  that  it 
is  actually  the  limiting  form  of  the  Lindman-Drummond  (Ref.  13)  ion  acoustic 
beam-whistler  type  instability. 


colllsionless  perpendicular  (Ref.  10)  and  parallel  (Ref.  2)  shock  waves 
formed  ahead  of  expanding  high  altitude  nuclear  debris  pistons.  The  fol- 
lowing two  instabilities  have  been  found  to  be  especially  attractive  as 
collisionless  mechanisms  for  pickup  and  heating  of  air  electrons: 

(i)  Current-driven  modified  two-stream  instability  of  McBride  et  al 
(Refs.  7 to  9)  which  should  be  operative  over  most  of  the  debris 
bubble  surface  excluding  the  loss  cone  exits. 

(ii)  Ion  cyclotron  beam  mode-whistler  mode  (current-free)  instability 
of  Golden  et  al  (Refs.  1,  2)  found  to  be  operative  along  the 
field  lines  and  particle  trajectories  issuing  directly  from  the 
loss  cone  exits. 

We  have  modelled  the  collisionless  shock  waves  as  Mott-Smith  layers 
of  interpenetrating  unshocked  (background  air)  and  shocked  (thermalized 
debris-air  piston)  flows,  so  that  these  layers  are  natural  environments 
for  streaming  instabilities.  Our  investigations  (based  on  linear  dispersion 
theory,  see  Ref.  10)  reveal  that  the  magnetic  counterstreaming  ion-ion 
instability  of  Papadopoulos  et  al  (Ref.  5)  cannot  structure  collisionless 
perpendicular  shocks  (formed  ahead  of  the  portion  of  the  bubble  which 
excludes  the  loss  cone  exits)  since,  in  the  reference  frame  of  the  shock 
front,  the  ensuing  magnetosonic  modes  are  blown  downstream  out  of  the  shock 
layer  before  they  can  grow  to  sufficiently  large  amplitude  to  scatter  in- 
coming air  particles . It  is  far  more  likely  that  perpendicular  shock  layers 
are  structured  by  the  turbulence  arising  from  the  current-driven  modified 
two-stream  instability  of  McBride  et  al.  A published  reprint  of  our  recent 
work  on  the  crossfield  magnetosonic  two-stream  instability  and  its  role  in 
perpendicular  shock  wave  structures  (Ref.  10)  is  also  included  in  Appendix  A. 
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Our  recent  linear  dispersion  theoretic  investigations  of  parallel 
shock  layers  (Ref.  2)  which  form  ahead  of  the  debris  plasma  issuing  from 
the  loss  cone  exits  reveal  that,  for  a given  Mach  number  M > M*  » 2.77 

A 

(M  ■ V /C  , V = shock  velocity,  C ■ Alfvdn  speed  in  unshocked  plasma), 

A U n U A 

the  shocked  hot  ion  cyclotron  beam  mode  can  always  drive  unstable  a par- 
ticular whistler  mode  which,  in  the  rest  frame  of  the  shock  front,  is 
stationary  near  the  leading  edge.  An  analysis  of  the  shock  interior  re- 
veals that  the  shock  Mach  number  determines  the  portion  of  the  shock 
thickness  in  which  unstable  whistlers  are  stationary  in  the  shock  rest 
frame.  For  M = 2.77,  such  modes  may  stand  only  at  the  leading  edge, 

A 

whereas  for  stronger  shocks  (M  > 2.77),  they  may  stand  at  all  points 
between  the  leading  edge  and  some  interior  point  which  is  dependent  on 
shock  strength.  For  very  strong  shocks  (M  >>  1),  fully  one-third  of  the 

A 

shock  thickness  is  filled  with  these  modes,  which  can  therefore  grow  to 
large  amplitude  and  couple  the  quiescent  background  air  plasma  to  the 
expanding  debris-air  piston.  In  an  earlier  study,  where  the  shocked  ion 
cyclotron  beam  mode  was  less  realistically  modelled  as  cold  fluid  (Ref.  1) , 
M*  was  found  to  equal  5.5  suggesting  that  the  turbulent  whistler  mode 
mechanism  for  momentum  coupling  was  operative  only  in  the  intermediate 
strength  parallel  shock  structures.  Our  new  theory  based  on  the  more 
realistic  hot  ion  cyclotron  beam  model,  however,  predicts  that  this 
coupling  mechanism  is  also  operative  in  the  weaker  parallel  shocks  en- 
countered in  high  altitude  nuclear  explosions.  A more  detailed  discussion 
of  this  analysis  is  presented  in  Sec.  B. 
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B.  THEORETICAL  ANALYSIS  OF  LOW  FREQUENCY  TWO  STREAM  INSTABILITIES 
AMD  THEIR  ROLE  IN  COLLISIONLBSS  SHOCK  STRUCTURE 


The  material  in  this  section  comprises  a manuscript  submitted  for 
publication.  It  summarizes  our  theoretical  work  on  low  frequency  two  stream 
instabilities  and  their  application  to  parallel  shock  wave  structure. 


■- 
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I.  INTRODUCTION 


It  has  been  theoretically  demonstrated  that  low  frequency  (w  1-5 
= ion  cyclotron  frequency)  instabilities  can  arise  from  relative  ion 
streaming  motions  along  a constant  applied  magnetic  field.1*2  These  were 
first  analyzed  by  Kovner1  who  considered,  among  other  current-free  initial 
flow  configurations,  one  where  (in  the  rest  frame  of  the  electrons)  a tenuous 
ion  beam  penetrates  a relatively  dense  ion  plasma;  here,  the  tenuous  ion 
cyclotron  beam  mode  interacts  with  well-defined  right  circularly  polarized 
ion  cyclotron  and  whistler  modes.  For  the  case  of  parallel  propagation 
(k  | | B^,  ^ = constant  applied  magnetic  field)  through  cold  plasma,  Kovner 
showed  that  only  the  whistler  modes  (corresponding  to  interactions  B and  C 
in  Fig.  1)  can  become  unstable.  More  recently.  Golden,  Linson,  and  Mani2 
extended  the  cold  plasma  part  of  Kovner's  work,  to  take  account  of  ion  beams 
having  densities  comparable  with  the  plasma  density. 

The  most  significant  result  reported  by  Golden  et  al,  however,  is  that 
these  unstable  ion  beam-whistler  wave  interactions  could  produce  the  turbulent 
structure  of  parallel  shock  waves  (shock  waves  which  propagate  along  a 
constant  applied  magnetic  field)  in  collisionless  plasma.  By  solving  the 
jc  | | linear  dispersion  relation  at  points  throughout  the  shock  layer 
under  the  assumption  that  the  local  ion  distribution  function  is  a Mott-Smith 
superposition  of  interpenetrating  unshocked  and  shocked  ion  flows  (electrons 
were  treated  simply  as  warm  fluid),  the  authors  found  that  certain  unstable 
modes  near  the  leading  edge  of  the  shock  could  propagate  at  precisely  the 
velocity  of  the  leading  edge  for  Alfven  Mach  numbers  MA  5.5.  This  suggests 
that  such  modes  have  ample  time  to  grow  to  sufficiently  large  amplitude  to 
scatter  incoming  (unshockcd)  ions  and  create  the  required  dissipation  for 
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Fig.  1 


Schematic  of  the  linear  dispersion  relation  for  parallel  propa- 
gation  showing  the  possible  interaction  of  the  right  circularly 
polarized  ion  cyclotron  and  whistler  branches  with  the  drifting 
ion  cyclotron  mode. 
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intermediate  strength  shocks.  Their  theory  has,  however,  two  serious  limi- 
tations imposed  for  the  sake  of  mathematical  simplicity:  (i)  only  unstable 

modes  propagating  along  B were  considered  and  (ii)  the  interpenetrating  ion 
flows  were  unrealistically  modelled  as  cold  raonoenergetic  beams.  Concerning 
the  first  of  these,  if  one  considers  waves  propagating  at  an  angle  to  B , 
electron  thermal  effects  have  a significant  interaction  with  the  magnetic 
modes.  A discussion  and  analysis  of  the  resulting  dispersion  relation, 
which  is  much  more  complicated,  is  deferred  to  a later  paper.  Concerning 
the  second  limitation,  it  is  at  once  apparent  that  the  shocked  ions  are,  by 
definition,  therraalized . Moreover,  even  the  unshocked  ions  must  be  suffi- 
ciently warm  (C,,  = -\J 5tcT/3 > C ^ , n^  = unshocked  ion  density) 

to  preclude  the  possibility  of  "switch-on"  type  shocks3  (rotation  of  the 
magnetic  field  across  the  shock  layers). 

The  main  objective  of  the  present  paper  is  to  re-examine  the  role  of 
these  ion  streaming  instabilities  in  parallel  shock  wave  structures  by  more 
realistically  taking  account  of  the  ion  thermal  effects.  This  we  do  first 
by  allowing  the  shocked  ions  to  be  Maxwellian  with  final  density,  mean  velocity, 
and  temperature  dictated  by  the  gas-dynamic  Rankine-Hugoniot  relations.  Then 
following  the  Mott-Smith  formalism,  the  ions,  in  the  rest  frame  of  the  shock 
front,  are  modelled  as  the  superposition  of  two  streams:  hot  downstream 

particles  drifting  through  the  cold  upstream  particles.  We  shall  see  that 
for  M > 2.77,  the  shocked  hot  ion  cyclotron  beam  mode  can  always  drive 

A 

unstable  a particular  whistler  mode  which,  in  the  rest  frame  of  the  shock 
front,  is  stationary  near  the  leading  edge.  Moreover,  such  stationary  growth 
modes  apparently  exist  up  to  one-third  of  the  distance  across  the  shock 
layer  starting  from  the  leading  edge.  Thus  the  turbulent  whistler  mode  momentum 
coupling  mechanism  for  intermediate  strength  shocks  (MA  >5.5)  suggested  by 


the  previous  cold  ion  cyclotron  beam  theory  is  suggested  as  well  by  the  hot 


beam  theory  of  the  present  paper  for  the  weaker  parallel  shocks  in  the  range 
2.77  < MA  < 5.5.  Finally  we  discuss  the  effects  of  relaxing  the  cold  up- 
stream plasma  assumption. 

Before  going  iito  the  shock  wave  problem,  it  is  instructive  to  first 
analyze  the  linear  dispersion  relation  in  the  initial  (current-free)  con- 
figuration  where  the  warm  electrons  have  zero  mean  velocity,  and  where  the 
hot  beam  ions  counterstream  through  the  cold  plasma  ions  along  the  constant 
external  magnetic  field.  This  we  do  in  Section  II.  The  application  of 
these  ion  streaming  instabilities  to  parallel  shock  wave  structure  is  then 
carried  out  in  Section  III  where  now,  in  the  rest  frame  of  the  shock  front, 
the  mean  motion  of  the  warm  electron  fluid  in  the  shock  layer  is  governed  by 
the  requirements  of  local  charge  and  current  neutrality.  Finally  in  Section 
IV,  conclusions  are  drawn  and  recommendations  arc  made. 


II.  ELECTRON  FRAME  ANALYSIS 

Consider  an  unbounded  collisionless  plasma  consisting  of  two  ion  beams 
of  densities  nj  and  n2  counterstreaming  with  velocities  Vj  and  V2  parallel 
to  a steady  magnetic  field  taken  to  lie  along  the  z axis.  We  assume  the 
unperturbed  background  electrons  to  be  described  by  a zero  mean  Maxwellian 
distribution  function.  The  charge  neutrality  and  zero  current  conditions 
are  given  by 

nj  + n2  - n (1) 

e 


and 


nlJ£l  + n2^2 


0 


(2) 


where  n^  is  the  number  density  of  electrons.  For  the  low  frequency  case  that 

we  are  interested  in,  it  is  crucial  that  the  zero  current  condition  (2)  be 

met  for  arbitrary  values  of  the  beam  strength  (0  < n2/n  « 1).  Deviations 

e 

from  the  unperturbed  state  are  assumed  to  be  governed  by  the  linear  Vlasov 
equation  with  an  unperturbed  ion  distribution  function  given  by 

fi<X>  = nl<5(v-^i)  + n2(irC2)-3/2  exp  [- (v-^>) 2/c2j  (3) 

where  c|  = 2<T2/m^. 

For  the  case  of  parallel  propagation  (jc  | | J!q)  the  linear  dispersion 
relation  is,  in  the  electron  rest  frame, 

* %,<*•“>  ■ <^>2-  <*> 


The  dielectric  components  in  (4)  are  calculated  by  addition  of  the  ion  and 

electron  polarizabilities.  These  well  known  polarizabilities  are  formulated 

explicitly  in  terms  of  the  (unperturbed)  ion  and  electron  distribution  functio  is 

and  therefore  take  account  of  the  bimodal  nature  of  the  ions  [Eq.  (3)]  and 

the  locally  Maxwellian  behavior  of  the  electrons.  We  make  the  following  low 

frequency  approximations:  (i)  the  vacuum  displacement  current  is  negligibly 

small  (u  <<  kc) , and  (ii)  the  electrons  have  zero  mass  (liol  <<  fi  = eB  /m  c) 

ce  o e 

so  they  can  be  modelled  as  a warm  fluid  (kr^e  <<  1).  The  dispersion  relation 
in  the  electron  rest  frame  ultimately  becomes 


D(k,to)  ■=  (1-n) 


ci 


w-kVi+ft^  Vrrc2 


r <ficl"ku 

fire  2 / (to- 

' — 00 


-ku)  exp[-(u/c2)2]  du 


kv2+ncl)-ku 


k2C2  u 

-i  - — f + — 

fic?  flci 


(5) 
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* 


Here  n “ n.,/ne  a measurti  °f  beam  strength.  Note  further  that  k,  V]  , and  V2 
art  the  respective  z components  of  V[,  and  V2  (and  are  tlius  signed  quantities). 
It  will  also  be  useful  to  incorporate  the  relati"e  drift  velocity  U into  (5) 


where 


U 


v2  - Vi  , 


-nU  , v2  = (l-ri)u 


(6) 


and  Eqs.  (1)  and  (2)  have  been  used.  Note  that  in  the  limit  c2  -*•  0 we  recover 
the  cold  dispersion  relation  of  Ref.  2.  The  solution  of  Eq.  (5)  w *=>  w(k) 
will  provide  the  spectrum  of  unstable  growth  modes  for  current  and  charge 
neutral  streaming  configurations  parametrized  by  n and  U.  We  separately 
consider  the  cases  of  small  and  arbitrary  beam  strengths  (i.e.  n <<  1 or 
0 < q < 1) , with  the  weak  beam  analysis  providing,  a simple  check  of  the  arbitrary 
beam  strength  calculation,  os  well  as  an  explicit  expression  for  the  linear 
growth  rate. 

A.  Weak  Beam  Analysis  (n  <<  1) 


In  the  weak  beam  approximation  we  set  u>  = v + iy,  y <<  |vl  and  perform  a 
perturbation  expansion  of  Eq.  (5)  for  y and  n small.  We  find  that  to  zeroth 
order  in  n there  is  no  contribution  to  y but  the  real  part  of  D(w,k)  = 0 leads 
to  the  well  known  whistler  dispersion  relation 

1/2 


v (k) 
o 


k2C2 


n 


ci 


2“ci 


kC 


A + A 


ci 


1 k2C2 


1 + - 


4 flc| 


(7) 


The  two  branches  of  this  dispersion  relation,  VQf,  are  shown  schematically 
in  Fig.  1.  To  obtain  the  first  order  contribution  to  the  growth  rate,  we  first 
rewrite  the  integral  term  of  Eq.  (5)  as 
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n e_x  (slcl/kc2  - x) 

- — dx 

ViT J (x-xq)  - iy/kc, 


x * u/c^ 

xo  a (v0-kU+Qcl)/kc2  . 

For  the  initial  value  problem  to  be  well-defined,  the  behavior  of  the 
integral  term  in  Eq.  (8)  for  y = 0 is  interpreted  as  the  limit  in  which 
y -»  0 through  positive  values.  For  y < 0 the  analytic  continuation  is  used  by 
deforming  the  path  of  integration  into  the  Landau  contour5.  To  account  for 
possible  changes  in  the  sign  of  k,  we  evaluate  (8)  using  the  Plemelj  formula 


f (x) dx 
(x-x0)  - ie 


P ~~~  + ffi(sgnk)f  (xQ) 

J-OD  O 


(where  c = y/kC, , y > 0)  and  obtain  to  first  order  the  growth  formula 


y VTrn  (vo+^ci^2  (kU-  \>Q) 

«ci  I k I c2  v0(vo+2nci> 


exp  - 


(v0-kU+flci)2 


where  we  have  used  |kj  = k(sgnk). 

Since  our  objective  in  analyzing  the  electron  frame  results  for  this 
instability  is  eventually  to  model  a collisionless  shock  wave,  we  consider 
streaming  configurations  for  which  U < 0.  It  is  clear  that  strong  interactions 
between  the  weak  warm  beam  and  the  cold  plasma  occur  in  the  vicinity  of 
intersections  between  the  "beam"  line 


vb  - kU  - ncl 


und  the  whistler  dispersion  vQ.  Since  U is  negative,  the  intersections  near 
which  growth  may  occur  are  shown  in  Fig.  1 to  be  with  the  negative  whistler 
branch  (i.e.  k < 0,  vQ_  > 0).  From  (10)  we  obtain  a condition  for  marginal 
stability 

Y = 0 when  vQ  = kU.  (11) 

Instability,  moreover,  can  occur  only  for  kU  > vD.  These  results  are  discussed 
in  the  context  of  arbitrary  beam  strengths  in  the.  next  section. 

B.  Arbitrary  Beam  Analysis 

In  this  section,  we  consider  beam-plasma  interactions  for  the  general 
case  when  the  strength  of  the  hot  beam  is  arbitrary  (i.e.  n is  no  longer 
small).  This  is  tantamount  to  an  investigation  of  growth  modes  in  the  in- 
terior of  a collisionless  shock  wave  modelled  as  two  interpenetrating  ion 
beams  with  beam  strength  n parametrizing  position  in  the  shock  wave  (Sec.  III). 
The  weak  beam-plasma  interaction  serves  to  model  the  leading  edge  of  such  a 
shock  wave. 

We  begin  our  analysis  by  observing  that  for  a given  beam  strength  n> 
the  stability- instability  boundary  is  represented  by  a curve  in  the  U-k  space. 
We  derive  a family  of  such  curves  by  finding  the  region  of  U-k  space  in  which 
the  solution  to  the  dispersion  relation  (5),  w(k),  has  a positive  imaginary 
part,  using  the  methods  of  complex  analysis.  In  particular,  the  well  known 


Nyquist  method  (5)  is  based  on  noting 


where  N - P represents  the  difference  between  the  number  of  zeros  and  poles 
of  the  complex  function  D(w)  (each  counted  according  to  multiplicity)  in  the 
region  of  the  w-planc  bounded  by  the  closed  contour  O.  The  right-hand  side 
of  the  equation  represents  the  change  in  the  argument  of  D(w)  as  C is  traversed 
once  in  the  counterclockwise  direction  in  the  w-plane.  We  note  that  D(w)  is 
analytic  off  the  real  axis,  so  that  if  one  chooses  C to  be  a semicircle  in 
the  upper-half  co-plane  (Fig.  2a) , then  one  may  set  P = 0.  Thus  we  may  write 


L£U 

+ 2 it  XI 


where  the  contribution  to  [Arg  D]  has  been  separated  into  the  two  parts  from 
Cj  (the  semicircular  arc  on  which  to  = Re-*-®)  and  '?jj  (the  real  axis  for  which 
-R  < w < R) , as  shawn  in  Fig.  2a.  In  the  limit  as  R + ”,  C encloses  the  entire 
upper-half  w-plane,  so  that  the  existence  of  unstable  whistlers  will  be  pre- 
dicted by  the  value  of  N,  the  number  of  roots  of  D = 0 with  positive  imaginary 
part.  VJe  introduce  here  the  analytic  function  Z(X/k)  defined  for  ImX  > 0 as 


Z(X/k)  ■=  it"1/2  / e— - 

I x - X/k 


and  as  its  analytic  continuation  for  ImX  s<  0 as  prescribed  by  the  Landau 
contour  . Evaluation  of  this  function  for  ImX  » 0 (X  = Xr  + iX^;  Xj^  = 0)  is 
obtained  by  using  the  Plemelj  formula  Eq.  (9)  and  yields 


r 00  _x2 

Xr/k)  ■ it-1/2  Pf  — - + u1^2i(sgnk)e 

J_oo  x-Xr/k 


-X2/kz 


Now  observe  that  the  integral  in  Eq.  (5)  may  be  written  as 
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-PLANE  D-PLANE 


-(u/ c z) 


2 


tict-kvi 

(u-kV2+fici) 


1 + 


u-kV2 


kc- 


w-kV2+S2ci' 
kco 


Now  for  R -*■“>,  wc  note  from  properties  of  7.  that 


D(w)  u “ ^ e^®  , 

fici  ^ci 


0 < 0 < tt  ; 


thus  [Arg  D]Cj  B it.  Therefore,  in  order  that  N be  nonzero,  we  must  require 
that 


[Arg  D]C 


II 


n(2n  + 1)  ; n * 0,1. . , 


First  we  split  D into  its  real  and  imaginary  parts  (for  u>  real): 


ReD  = (1-n) 


«ci 

oi-kV2^ 

Rezj 

f vjJ-kV2+nci ) 

1 

k2Cj[  u> 

- + n 

u- kV  i -Hi  . 

1 Cl 

IK 

CM 

U 

_1  ’ ne|  + oci 

ImD  = 


o-kV2 

kc2 


IraZ 


-kv2+Qci\ 


kc; 


We  now  note  that  for  w -*■  ±®  we  have 

ReD  -v  to/nci  -*■  ±“ 


and 

ImD  'v.  tt1/2  (w/ 1 k j C2)  exp(-co2/k2c2)  •+  ±e 

so  that  the  point  representing  D in  the  complex  D-plane  traverses  a curve 
from  - » - it  ->•  + °°  + ie  as  w traverses  O^.  From  Eq.  (17)  we  see  that  this 
point  crosses  the  real  D axis  once  only,  at  « ■ kV2  so  that  three  distinct 
curves  in  the  D-plane  may  be  traced  as  u traverses  Cjj*  These  are  shown  in 


(16) 


(17) 
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Fig.  2b  as  A,  B,  and  C and  are  keyed  to  the  sign  of  ReD  at  the  point  for 
which  ImD  = 0.  In  summary,  for  u>  m kV2  and  ImD  ■ 0, 

A ReD  > 0 =>  [Arg  DJC^^  = = 1 (Instability), 

B ReD  < 0=>[Arg  D)c  - -*=>N  “ 0 (Stability), 

C ReD  “ 0 "=^>  (Marginal  Stability). 


Thus  we  see  that  the  domain  of  instability  will  be  determined  by  the  condition 
that  ReD  >,  0 for  w » kV2 , with  the  equality  sign  giving  the  stability-insta- 
bility boundary  in  the  U-k  plane  for  an  arbitrary  beam  strength  n.  Now 
substituting  u * kV  into  ReD  gives  the  instability  criterion 


(l-n) 


k2u2/nc| 

1 - k|uj/Dc. 


k2C2 


ci 


> 0 


(lb) 


where  |u|  =■  -U  for  the  negative  drift  case  considered  here.  Two  observations 
are  immediately  seen:  (i)  for  k|u|  > nci  (k  > 0),  there  will  be  unconditional 

stability  for  all  beam  strengths,  and  (ii)  for  n -*•  1,  there  can  be  no  insta- 
bility unless  U -*•  ® and  k < 0.  Thus  for  a cold  beam-hot  plasma  configuration, 
instability  will  only  occur  for  modes  propagating  antiparallel  to  BQ  generated 
by  a large  relative  streaming  velocity.  Observations  (i)  and  (ii)  thus  iraply 
that  unstable  whistlers  may  exist  for  finite  U only  if  |u|  < 0ci/k  and  n < 1. 
With  these  constraints  on  |u|  and  n the  relevant  inequality  for  instability 
becomes 


Mp  > H(x,n) 


where 


H(x,n) 


2 (l-n) 


l-n  4 (l-n) 2 


1/2 


(19) 


(20) 
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E 


Md  - |u|/cA  , 


X 


“ kCA^ci  * 


A simple  asymptotic  analysis  of  Eq.  (20)  shows  that 


lira  H(x,n) 
n-+l 


l/x 


x > 0 


x < 0 


so  that  as  q varies  between  0 and  1 the  curves  of  marginal  stability, 

MD  = H(x,q) , fill  the  region  between  = H(x,0)  and  = l/x  (x  >0).  Thus 
for  a given  beam  strength,  the  domain  of  MQ-x  pairs  for  which  instability  may 
occur  is  given  by  the  inequality 


H(x,n0)  ^ ^ l/x  • 


and  is  represented  graphically  as  the  region  between  the  curves  H(x,n0)  and 
l/x.  Plots  of  H(x, q) , for  various  values  of  q,  end  l/x  are  shown  in  Fig.  3. 
Note  that  for  q <<  1 the  marginal  stability  curve  is  simply  related  to  the 
negative  branch  of  the  whistler  dispersion  curve  vQ_  by 


M 


D 


1/2 


1 Vq- 
x Qci 


in  agreement  with  the  results  of  Sec.  II. A.  and  F.q.  (11). 


III.  APPLICATION  TO  PARALLEL  SHOCK  WAVE  STRUCTURE 


The  analysis  of  the  preceding  section  can  now  be  applied  to  the  structure 
of  parallel  shocks.  We  determine  both  the  spectrum  of  unstable  modes  that  may 
exist  at  various  points  in  the  interior  of  the  shock  as  well  as  which  of  these 
modes  are  stationary  in  the  rest  frame  of  the  shock  leading  edge.  In  this 
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work  the  cold  ion  stream  requirement  of  the  previous  theory2  is  partially 
relaxed  by  allowing  the  shocked,  downstream  ions  to  have  a Maxwellian  velocity 
distribution  with  temperature  and  mass  velocity  dictated  by  the  conser- 
vation relations  applied  across  the  shock  wave.  The  unshocked,  upstream  ions 
are  assumed  to  be  very  cold  compared  to  the  downstream  ions  and  are  moving 
with  mass  velocity  V^.  The  shock  layer  is  now  modelled  as  a superposition 
of  these  up-  and  downstream  ion  beams  moving  with  relative  velocity  JJ  = ^ - Vy 
and  with  beam  streng-.h  n(z)  dependent  on  position  in  the  shock  layer.  Using 
the  Mott-Smith  assumption  that  the  velocity  distribution  in  the  interior  of 
a shock  wave  be  a bimodal  Maxwellian  then  leads  to  the  ion  distribution  in 
the  shock  layer  in  tae  form 

fi(z*^  = au(z)5W+nd(z)(,,cd)'3/2  exp[-(ry2/cd] 

where  C2  = 2<T ,/ni . (<  = Boltzmann's  constant)  and  n ,(z)  represent  the  up- 
d d 1 u,d  r 

and  downstream  ion  densities  respectively  at  a point  z in  the  shock  layer. 

The  warm  background  electrons  move  with  mass  velocity  ^g(z)  = [l-n(z)]^u  + n(z)^ 
where  n(z)  = nd(z)/[nu(z)  + n^z)]  dictated  by  the  requirements  of  local  charge 
and  current  neutrality.  For  simplicity  n^(z)  is  assumed  to  decrease  linearly 
front  its  upstream  value  n at  the  shock  leading  edge  to  zero  at  the  trailing 
edge  whereas  n^z)  is  assumed  to  increase  linearly  from  zero  at  the  shock 
leading  edge  to  its  downstream  value  n^  at  the  trailing  edge.  A schematic 
of  the  density  distribution  is  sketched  in  Fig.  4?. 

The  valid  dispersion  relation  is  now  obtained  from  Eq.  (5)  by  Doppler 
shifting  the  frequency  to  according  to 

to  = to  + k-V  ( 21a) 

s ‘V  "Ve 

and  by  transforming  all  velocities  from  the  electron-  to  shock-frame  according 
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(HOT  SHOCKED  PLASMA) 


to 


V„  “ v,  + V„ 

*\,u  i *v,e 


V , = V + V 

'vd  'v 2 've 


By  substituting  those  into  Eq.  (5)  we  obtain 


(21b) 

(21c) 


D(w  ,k)  = [l-n(z)] 


‘ci 


-1  - 


3-kVu+«ci 

k2ci(z) 


- + n(z) 


1.  + 


Vs-kVdl 

kC,: 


/“  s-kVd 

Wci| 

kCd 

\ 

/ 

*ci 


Ws  kVft(z) 
^ci  ^ci 


(22) 


The  solution  of  Eq.  (22),  ws  = tos(k),  will  provide  the  spectrum  of  unstable 
growth  modes  for  a particular  shock  position  z and  given  shock  strength, 

Vu  - Vd.  Only  those  growth  modes  which  are  stationary  in  the  shock  rest 
frame  (i.e.  have  zero  group  velocity)  are  of  interest,  since  only  these  modes 
have  sufficient  time  to  grow  to  large  amplitude  to  scatter  incoming  ions. 

The  analysis  of  growth  and  group  velocity  of  these  modes  is  again  separated 
into  two  parts:  (i)  the  shock  leading  edge  (p  <<  1),  and  (ii)  arbitrary 

points  in  the  shock  interior  (0  < p < 1) . In  a«dition,  the  leading  edge  calcu- 
lation has  been  extended  to  include  the  effects  of  finite  temperature  in  the 
unshocked,  upstream  ions. 

A.  Leading  Edge  Analysis  (p(z)  <<  1,  Cold  Upstream  Plasma) 

The  leading  edge  analysis  proceeds  in  the  same  way  as  the  weak  beam 
analysis  of  Sec.  II.  Setting  u>g  = vs  + iys,  0 < ys  <<  |vs|  and  performing 
a perturbation  expansion  for  both  ys  and  p small  leads  to  the  solution 
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(23a) 


v8(k;V  " kVu  + k2CA2/2«ci  - kCAo[l  + (kCAo/2nci)Z)1/z  , 


Y(k;V  ) 

(v-kVu+fic1)2(kVd-v) 

(v-kVj+ficl)2 

nci  lilt'd 

(v-kVu) (v-kVu+2ncl) 

k2C2 

(23b) 


Strong  interaction  evidently  corresponds  to  the  intersection  of  the 
v - kVj  + S2C^  = 0 hot  ion  beam  line  with  the  cold  whistler  dispersion  curve, 
Eq.  (23a),  in  the  third  quadrant  of  the  (ws,k) -plane.  Here  the  resulting 
Ymax  varies  directly  with  the  strength  q of  the  hot  beam.  In  the  previous 
theory2,  where  the  shocked  ions  were  modelled  as  a cold  energetic  beam,  the 
corresponding  weak  ion  beam-whistler  interaction  was  found  to  be  stronger, 
i.e.  YmaX“Vn . Thi*  is  not  surprising  since,  in  the  present  theory,  the 
whistler  wave  interacts  strongly  only  with  a small  percentage  of  the  hot 
beam  ions,  whereas  in  the  previous  theory,  all  of  the  particles  in  the  mono- 
energetic  ion  beam  participate  in  the  interaction. 

From  Eq.  (23b),  we  see  that  instability  (i.e.,  y > 0)  is  possible  only 
if  kV^  - v(k;Vu)  > 0;  this  together  with  (22)  gives  y > 0 for 


Nca 

2«ci 


k2<’  2 
, k l'Ao 
+ y— 

AO  ? 

ql‘Ci 


1/2 


(24) 


V^,  however,  is  connected  to  Vu  through  the  gasdynamic  Rankine-Hugoniot 
relation 


Vd  M2  + 3 

— - • (HA  “ VCAo>  Cp/Cv  * 5/3)  <25) 

vu 


where,  to  avoid  the  possibility  of  "switch-on"  type  shocks3,  we  have  assumed 


CAo“  Cu  = 5xTu/3m£.  Equations  (24)  and  (25)  combine  to  give  the  leading  edge 


instability  condition 


M*  > 


|g  + (1  + |g2) , 


(26) 


8 = -S  + (l+^.)l/2  , x - ^ , CA0  - CA(0)  . 


Q 


ci 


Next  from  Eq.  (23a)  the  zero  group  velocity  condition  is  easily  found  to  be 


M,  = -x  + (1  + *i)1/2  + *i(l  + *1) -1/2  . (27) 

a 4 4 4 

A sketch  of  the  stability-instability  boundary  [Lq.  (26)]  and  the  locus  of 
zero  group  velocity  inodes  [Eq.  (27)]  is  shown  in  Fig.  5 plotted  in  the 
(Ma,  | x ! ) plane.  Evidently,  for  a given  Mach  number  MA  >2.77,  the  shocked 
hot  ion-cyclotron  beam  mode  can  drive  unstable  a particular  whistler  mode 
which,  in  the  rest  frame  of  the  shock  front,  is  stationary  near  the  leading 
edge.  This  whistler  can  therefore  grow  to  large  amplitude  there  and  ulti- 
mately couple  the  upstream  and  downstream  plasma  flows.  Thus  the  turbulent 
whistler  mode  strut ture  for  intermediate  strength  shocks  (MA  >5.5)  suggested 
by  the  previous  cold  beam2  theory  is  suggested  as  well  for  weak  shocks 
(2.77 < Ma  < 5.5)  by  the  present  hot  beam  theory. 

B.  Effects  of  Finite  Upstream  Temperature  on  the  Structure  of  the 

Shock  Leading  Edge 


These  leading  edge  conclusions  are  modified  somewhat  by  including  in 
the  analysis  the  effects  of  finite  upstream  temperatures.  The  valid  disper- 
sion relation  in  this  case  contains  two  singular  integrals  of  the  form  found 
in  Eq.  (22).  To  obtain  some  indication  of  temperature  effects  in  the  unshocked 
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plasma  thermal  speed,  Cu,  is  finite  but  sufficiently  small  so  that 

k2c2  « 2<vs  - kVu  + fici)2  . 


(28) 


This  approximation  results  in  limiting  lowest  order  upstream  thermal  effects 
to  the  imaginary  part  of  D(ws,  k)  only.  Thus  the  zeroth  order  approximation 
to  es  is  again  given  by  the  whistler  dispersion  formula,  Eq.  (23a).  In 
addition,  it  is  assumed  that  in  the  imaginary  part  of  D(ws,  k) , both  damping 
and  growth  effects  occur  to  the  same  order,  leading  to  a growth  rate  y of 
the  form 


•ci 


, : ^s+|k|vd+ficl)?- 

, (n-nc)  -j^j^r  [''s+IMVdl  exp  j - 


k2cg 


(29) 


and  with  the  instability  criterion  (k  < 0) 


vs  < ~ lkl Vd* 


(30) 


which  is  identical  to  that  given  by  Eq.  (24).  Consequently,  the  stability- 
instability  bound;-  ry  is  again  given  by  Eq.  (26)  is  shown  in  Fig.  5.  Further- 
more, in  Eq.  (29),  nc  represents  a critical  penetration  distance  into  the 
shock  before  which  (n  < nc)  the  effects  of  damping  dominate  growth  but  beyond 
which  instability  again  occurs.  Some  simple  calculations  indicate  that  for 
Ma  > 3,  Eq.  (28)  is  well  satisfied  and  nc  < .05,  so  that  the  approximate  cal- 
culation is  consistent  with  the  leading  edge  small  perturbation  theory. 
Consequently,  the  ’ffects  of  finite  upstream  heating  (sufficient  e.g.  to 
prevent  the  occurrence  of  switch-on  shocks,  but  small  enough  so  that  Eq.  (28) 
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is  sa  ed)  is  to  cause  a small  portion  of  the  shock  interior  to  be 

stable  to  the  two  stream  disturbance.  The  conditions  for  the  instability 
to  occur  are  then  both  Eq.  (30)  and  n > 17.  For  the  range  of  shock  para- 
meters of  interest,  however,  nc  < -05  so  that  the  principle  results  of  this 
section  are  unaffected. 

C.  Shock  Layer  Analysis  (0  < n < 1,  Cold  Upstream  Plasma) 

The  previous  analysis  of  the  shock  leading  edge  has  demonstrated 
that  relatively  weak  shocks  (M  z 2.77)  can  be  structured  by  this  Whistler 

A 

mode  instability  since  for  all  shock  Mach  numbers  larger  than  2.77,  there 
can  exist  unstable  whistlers  that  are  stationary  (have  zero  group  velocity) 
in  the  shock  rest  frame.  The  results  of  Sec.  II  indicate  that  for  all  beam 
strengths  unstable  whistlers  can  be  generated  by  sufficiently  large  relative 
beam  drift.  Consequently,  it  can  be  inferred  that  at  all  points  of  the 
shock  interior,  unstable  whistler  modes  may  exist.  Since  the  instability 
criterion  as  given  in  Eq.  (18)  is  invariant  with  respect  to  the  transfor- 
mation of  Eqs.  (21),  the  criterion  for  the  shock  interior  may  be.  written 


k2U2/fl  2 

[1  - nOO) 7—7- 

1 - kU/B 


k2C2(z) 


where  luj  * V - V,  and  the  variation  of  n and  C,  with  position  in  the  9hock 
1 1 u d A 

layer  is  explicitly  noted.  Furthermore,  applying  Eqs.  (21)  to  the  Nyquist 
criterion  u>  = kV  leads  to 


w “ kV  , . 
s d 
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Therefore  for  points  on  the  boundary  of  instability  (equality  in  Eq.  (31)), 

<*>s  is  real  and  is  given  by  Eq.  (32)  for  all  positions  in  the  shock  layer. 
Again  it  may  be  observed  for  Eq.  (31)  that  instabilities  may  exist  for  finite 
drift  only  if  k|ll|  < (k  > 0)  and  n < 1 . Thus  trailing  edge  instabilities 
may  exist  only  for  infinitely  strong  shocks.  With  these  restrictions  on 
| U | and  n,  the  relevant  condition  can  be  recast  as 


r 1 i 

kC £<*>  , 

CA(Z) 

i 

K-* 

/-“V 

N 

2U 

ci 

l-n(z)  4[i-n(z)]2fic^ 

(33) 


with  the  equality  defining  the  surface  of  marginally  stable  modes. 

Now  within  the  framework  of  the  Mott-Smith  Formalism,  the  shock  is 
considered  to  be  a superposition  of  upstream  and  downstream  plasmas  with 
beam  densities  "^(0  and  n^(z)  at  any  point  in  the  shock  with  relative  drift 
velocity  | U | = > 0.  For  simplicity  the  density  in  the  shock  is 

assumed  to  be  linear  in  z so  that  n^  z)  may  serve  to  define  the  fraction 
of  penetration  into  the  shock,  £,  as 


n (z)  - n 
e u“ 


u®- 


where  n and  n,  are  the  plasma  densities  in  the  uniform  regions  before 

and  after  the  shock  and  ng(z)  the  total  plasma  density  at  the  point  z in 

the  shock.  Thus  c = 0,  1 refer  to  leading  and  trailing  edge  conditions 

respectively.  Using  the  definition  of  n(z)  = n,  (z)/n  (z)  then  leads  to 

a e 


n 

n + t(1  - n) 


(34) 


where  t = n^/n^  = V u/V ^ is  the  shock  compression  ratio  expressed  in  terms  of 
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M by  Eq.  (25).  Therefore, 

A 


l - n(z)  = (1-0 


1 + 3 i. 


"Ill 

M2  + 3 
A 


(35) 


Furthermore,  the  variation  of  Alfv£n  speed  through  the  shock  can  be  expressed 
in  terms  of  leading  edge  conditions  by  the  relation 
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1 - n(z) 

1 - ; 


(36) 


Since  the  relative  velocity  |u|  = vu  “ Vjj  may  be  referred  to  the  upstream 
Alfvdn  Mach  number  M^  through  the  Rankine  Hugoniot  conditions  for  a parallel 
6hock  wave  as 


V - V, 
u d 


Ao 


3(M2  - ]) 
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(37) 


Eq.  (33)  can  be  written  solely  in  terms  of  M^,  k,  and  ; as  follows: 


MA  > 3«  + 


i + 4g2 


1/2 


h F (k,  O 
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where 
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(38b) 


In  addition,  the  region  of  unconditional  stability  (k|u|  < f<ci  for  k > 0) 
can  be  rewritten  as 


u 
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= G(k) 


(39) 


Plots  of  F(k,  c)  and  G(k)  are  given  in  Fig.  6.  / simple  asymptotic  analysis 

of  Eqs . (38)  shows  that  for  5 ->  1 we  have  the  twc  results 

lira  F(k,  4)  = G(k) ; k > 0 

4-1 

lim  F(k,  C)  *°  k < 0 , 

4-1 

so  that  the  family  of  curves  given  by  F(k,  0 fills  the  region  between 
4 = 0 and  the  boundary  of  unconditional  stability  given  by  G(k). 

The  result  is  chat  in  the  M -k  space  of  Fig.  6 at  a particular  shock 

A 

position  C,  the  domain  of  unstable  modes  must  lie  between  the  corresponding 
t,  curve  and  the  curve  of  unconditional  stability.  Furthermore,  as  £ -*■  1, 
this  domain  vanishes  rapidly.  Our  principle  result,  however,  is  that 
although  unstable  whistlers  may  exist  at  all  points  of  the  shock  layer, 
very  high  are  required  for  such  instabilities  near  the  shock  trailing 
edge.  Although  soma  of  these  modes  will  be  stationary  (have  vanishing 
group  velocity  in  t ne  shock  rest  frame),  it  seems  certain  that  the  weakest 
shocks  will  have  unstable  whistlers  with  zero  group  velocity  only  near  the 
leading  edge. 

This  latter  conclusion  can  be  strengthened  by  considering  the  group 
velocity  of  unstable  modes  at  each  6hock  position  as  a function  of  shock 
Mach  number  MA>  In  particular  on  the  marginal  stability  surface  (MSS), 
for  which 

Ma  = F(k,  4)  (40) 
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and  (*)  denotes  quantities  evaluated  on  the  MSS.  After  substituting 
Eqs.  (41)  and  (42)  into  the  real  and  imaginary  parts  of  Eq.  (22),  we  find 
(i)  that  the  zeroth  order  terms  in  the  6 ordering  vanish  identically  after 
using  Eq.  (40),  and  (ii)  to  within  0(6)  we  have 


k*(V  -V  *)fi  , 

. . _ «L  6. £i 

(k*u-nci)2 


k*(V  _V  *) 

n — s — Hi 

Q . 
ci 


In  these  expressions  Kt  and  Ij  are  related  to  the  real  and  imaginary  parts 
of  7.  for  real  argument  by  the  equations 


ft  . ft  . 

R 1 . — Z ?— 

|k*cdl  r\|k*Cdi 


where  Z c + iZ^.  In  principle,  Eqs.  (44)  and  (45)  can  be  solved  to  give 

the  group  velocity  of  any  marginally  stable  mode,  V * = V *(k*,  c).  In 

E 8 

particular  the  locus  of  stationary,  marginally  stable  modes  is  determined 

by  setting  V * = 0 in  Eqs.  (44)  and  (45)  and  solving  for  the  resulting 
8 

k*  - £ pairs.  Equations  (25)  and  (35-40)  can  be  used  to  express  n(z),  V , 
and  V in  terms  of  k*  and  £.  Furthermore,  the  arguments  of  Rj  and  Ij  can 
be  written  as 


°ci  CAo 

lk*lCAo  Cd 


'k*lCAo  " 6 Td 


where  we  have  used  C 2 » 5<T  /3m.  to  avoid  switch-on  shocks.  The  ratio  T /X, 

Ao  u l u d 

can  be  obtained  via  momentum  and  energy  conservat ion  across  the  shock  as 


T.  15M^  + 42M2  - 9 * 

d A A 


Using  these  relations  in  Eqs.  (44)  and  (45)  leads  ultimately  to  an  equation 
of  the  form 


x*  e C(x*,  c)  , 


where 
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Ci(x*.  rj  = 
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Equation  (46)  has  been  solved  numerically  for  x*(0  by  a convergent  iterative 
procedure.  The  results  are  shown  as  a locus  of  points  on  Fig.  6.  Further- 
more, Fig.  7 shows  a plot  of  vs.  £ for  which  marginally  stable  modes  are 
stationary  in  the  rest  frame  of  the  shock.  Figures  6 and  7 demonstrate  that 
and  k*  both  become  unbounded  as  £ -►  1/3.  An  asymptotic  expansion  of 
Eqs.  (47)  for  M >>  1,  x*  >>  1 leads  to  the  expression 

A 

C(x*,  £)  'V  A(;)x*  + B(c)/x* 


where 


and 
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B = (1  - £)  - C)3/2+  (|)3/2(1  - C)5/2  - 1 1 . 
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Fig.  7 Range  of  Alfv£n  Mach  numbers  for  which  unstable 
whistlers  are  stationary  in  the  shock  frame  vs. 
shock  position  £. 
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The  solution  of  Eq.  (46)  then  becomes 


where  the  unphysicaL  positive  solution  has  been  discarded.  Note  that 
Eq.  (48)  becomes  unbounded  as  ^ -*■  1/3  from  below  and  that  no  real  solution 
for  k*  exists  for  z 1/3.  This  is  consistent  with  our  numerical  solution 
of  Eq.  (46)  as  shown  in  Figs.  6 and  7.  The  existence  of  these  stationary 
unstable  whistler  modes  for  0 < £ < 1/3  significantly  extends  the  previous 
cold  theory2,  for  which  unstable  whistlers  were  stationary  only  for  0 H < .1. 
This  result  strengthens  the  previous  conclusion  of  Golden  et.  al.2  that  this 
two  stream  instability  is  an  important  source  of  turbulent  dissipation  re- 
quired to  structure  a collisionless  parallel  shock  wave. 


IV.  DISCUSSION  AND  CONCLUSIONS 


The  instability  generated  by  two  ion  beams  counterstreaming  with 
relative  velocity  U and  beam  strength  n parallel  to  an  external  magnetic 
field  and  in  the  presence  of  a charge  and  current  neutralizing  electron 
background  has  been  analyzed  using  linear  dispersion  theory.  The  previous 
cold  beam  theory  of  Golden  et.  al.2  has  been  extended  by  allowing  one  of 
the  ion  beams  to  have  a finite  temperature  in  order  to  model  a collisionless 
shock  wave  more  realistically.  By  using  the  Nyquist  method5,  a family  of 
curves  in  the  normalized  U vs.  k space  has  been  obtained,  which  represents 
the  locus  of  marginally  stable  whistlers  parametrized  by  the  beam  strength  n 
(0  n < 1) . The  instability  for  n <<  1 may  be  interpreted  as  a hot  ion 
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cyclotron  beam  mode  driving  unstable  a cold  plasma  whistler  mode.  The 
detailed  analysis  for  n <<  1 results  in  a growth  rate  y « n in  contrast 
to  the  previous  cold  beam-cold  plasma  theory  in  which  a typical  hydrody- 
namic result,  y « n1^2,  was  obtained. 

A collisionless  shock  wave  has  also  been  modelled  by  using  a Mott-Smith 
bimodal  Maxwellian  assumption  for  the  ion  velocity  distribution  in  the 
shock  rest  frame  with  constant  up-  and  downstream  temperatures  throughout 
the  shock  wave  but  with  linearly  varying  densities.  Each  point,  C,  in  the 
shock  then  serves  a3  an  appropriate  environment  for  the  two  stream  instability 
under  consideration,  where  here  the  beam  strength  n(C)  is  a measure  of 
position  in  the  shock.  We  have  obtained  the  family  of  curves  of  marginal 
stability  parametrized  by  shock  position  5 in  the  vs.  k space.  Superim- 
posed on  these  curves  in  the  locus  of  marginally  stable  modes  which  are 
stationary  in  the  shock  rest  frame.  It  has  been  found  that  for  a given  M^, 
stationary  unstable  whistlers  may  exist  at  all  points  in  the  shock  from  the 
leading  edge,  £ = 0,  to  the  interior  point,  £ = < 1/3  as  given  by  Fig.  7. 

These  modes  have  time  to  grow  to  large  amplitude  and  will  determine  the 
turbulent  spectrum  for  a given  shock  strength  (M^) . From  Figs.  5 and  6 
we  can  determine  for  a given  the  wave  numbers  of  the  leading  edge  and 

modes.  We  have  found  that  for  all  M,  > 2.77,  these  two  wave  numbers 
M A 

differ  by  approximately  2%.  This  suggests  that  for  a given  shock  strength, 
the  turbulent  spectrum  has  a very  narrow  band  width  and  is  located  near 
the  wave  number  of  the  unstable,  stationary  leading  edge  mode. 

The  principle  differences  between  the  previous  cold  theory2  and  the 
present  work  are: 

(i)  weaker  leading  edge  growth  rate  characteristic  of  Vlasov  plasmas. 
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(ii)  existence  of  unstable  stationary  whistlers  for  weaker  shocks 


(2.77 « M ) than  in  the  previous  theory  (5.5  < M ), 

A A 

(iii)  the  region  of  instability  extends  to  arbitrarily  large  shock 
strengths  but  was  bounded  above  in  the  previous  theory,  and 

(iv)  unstable  and  stationary  modes  may  exist  for  a larger  range  of 

shock  positions  (0  £ < 1/3)  than  previously  (0  < c < .1). 

The  effects  of  finite  upstream,  temperature , T^,  have  been  considered 

in  a simple  way  since  to  avoid  switch-on  shocks3  it  has  been  assumed  that 

5kT  /in.  = C 2.  We  have  found  that  for  the  range  of  shock  Mach  numbers 
' u i Ao 

under  consideration  the  results  of  our  cold  upstream  theory  remain  essentially 
valid  even  for  a finite  but  small  (Sec.  III.B.).  The  effects  of  obliquely 
propagating  modes  (J^  \\  J3q)  and  nonlinear  saturation  on  parallel  shock  struc- 
tures arc  currently  under  consideration  and  will  be  reported  in  a later  paper. 
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A principal  objective  of  this  research  is  to  investigate  the 
effects  of  high  altitude  nuclear  burrcs  on  energetic  charged  particle 
injection  into  the  earth's  magnetosphere.  The  work  reported  in  Sec.  B 
has  identified  a particular  two  stream  instability  that  can  serve  as  the 
necessary  dissipative  mechanism  to  structure  a collisionless  shock  wave. 

This  instability  occurs  in  the  loss  cone  exits  of  an  existing  nuclear 
debris  bubble  as  the  expanding  debris  from  a second  detonation  passes  over 
the  ionized  ambient  air  plasma.  The  importance  of  this  result  is  that 
since  any  finite  amplitude  wave  passing  over  a stationary  plasma  has  a 
tendency  to  put  the  originally  stationary  plasma  into  motion,  the  passage 
of  a debris  structured  collisionless  shock  wave  through  the  ionized  ambient 
air  in  the  loss  cone  corridor  will  serve  to  accelerate  these  ions  along 
the  distended  geomagnetic  lines  and  thus  generate  particle  fluxes  into 
the  magnetosphere. 

In  order  to  apply  the  results  of  our  theory  to  the  calculation  of 
these  fluxes,  we  consider  the  environment  generated  by  two  high  altitude 
Starfish-like  nuclear  bursts  (Fig.  1) . We  show  first  that  for  these  para- 
meters unstable  whistlers  traveling  with  the  shock  exist  and  have  sufficient 
time  to  grow  to  large  amplitude  while  passing  over  roughly  607.  of  the 
loss  cone  corridor.  Therefore,  this  emphasizes  the  relevance  of  the 
theory  presented  in  Sec.  B.  to  the  deposition  of  energetic  particles  into 
the  magnetosphere. 


Explicit  calculations  of  the  number  of  such  particles  are  performed 
using  three  different  approaches.  It  will  be  seen  that  all  three 
approaches  yield  the  same  order  of  magnitude  although  based  on  different 
physical  assumptions  about  the  shock  dynamics  after  shock  formation. 
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II.  GROWTH  RATES  FOR  STARFISH  PARAMETERS 


Consider  the  multiple  burst  environment  as  shown  In  Fig.  1.  If  a 
second  Starfish-llke  nuclear  burst  occurs  Inside  of  an  already  existing 
nuclear  debris  bubble  (of  extent  ^ 1000  km) , a fraction  of  the  second 
bomb's  debris  will  be  Injected  directly  into  the  loss  cone  exits  of  the 
parent  bubble. 


Figure  1 


The  number  of  such  debris  particles  which  enter  a loss  cone  exit  depends 
critically  on  the  mass  and  orientation  of  the  second  bomb  casing  just  prior 


to  the  burst.  In  the  loss  cone  corridor,  the  debris  plasma  drives  a shock 
wave.  This  shock  is  structured  by  unstable  whistler  waves  which  stand  at  Its 
leading  edge.  To  verify  that  these  whistler  modes  can  grow  to  sufficiently 
large  amplitude  to  scatter  incoming  air  ions  (as  viewed  in  the  reference 
frame  of  the  shock  front),  we  require  that  the  non-dimensional  parameter 
Y i/u  > 1,  where  l ^ 500  km  is  a characteristic  length  of 


travel  of  the  debris  particles  In  the  loss  cone  corridor.  During  day-  . 
light  burst  conditions  the  density  of  oxygen  Iona  is  n0+  ~105-106,  the 
larger  value  reflecting  maximum  sunspot  activity.  For  an  ambient  field 
strength  B ^ 0.3  gauss,  this  corresponds  to  an  Alfven  speed  CA  ~ 163  - 
516  km/sec  in  the  ambient  air  plasma. 

For  typical  initial  casing  velocities  of  from  500  km/sec  to 
2000  km/sec  the  initial  Alfvdn  Mach  numbers  are  from  1.96  to  12.22.  A 
summary  of  the  spectrum  of  stationary  unstable  leading  edge  whistlers, 
calculated  using  Eq.  (27)  of  Sec  B.  is  given  in  Fig.  2.  It  is  seen 
there  that  a broad  range  of  wave  numbers  can  grow  to  large  amplitude 
during  daytime  (n’^lO^  cm“^)  bursts.  To  investigate  whether  these  modes 
can  achieve  these  large  amplitudes  in  the  loss  cone  exit,  we  must  cal- 
culate the  growth  rate  of  these  unstable  waves.  For  simplicity  we 
choose  Jlc  * 2 for  which  “ 4.12,  indicating  a realistic 

initial  debris  velocity  of  674  km/sec.  In  this  case  we  calculate  the 
linear  growth  rate  for  leading  edge  modes  using  Eqs.  (23)  of  Sec.  B and 
the  Rankine-Hugoniot  shock  wave  relations  to  be  Y - 167  n.  At  n - .2 
(shock  leading  edge)  we  find  y sju  m 10,  indicating  10  e-folding  times 
pass  as  the  shock  traverses  500  km  of  the  loss  cone  corridor.  Thus  the 
[k  | ■ 2 Dc1/Ca  unstable  whistler  modes  can  grow  to  sufficiently  large 
amplitude  after  traveling  500  km  along  the  loss  cone  corridor  to  pick 
up  background  air  plasma.  The  air  electrons  are  rapidly  accelerated  by 
ensuing  ion-electron  electrostatic  instabilities.  Since  the  calculation 
is  based  on  the  lr*t  r bound  of  MA,  it  is  clear  that  10  e-folding  times 
in  500  km  is  a conservative  estimate  since  stnnger  shocks  will  have 
larger  growth  rates. 
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III.  CALCULATIONS  OF  PARTICLE  DEPOSITION 


Perhaps  the  simplest  method  for  estimating  the  additional  particle 
deposition  due  to  collisionless  mechanisms  is  to  assume  that  the  loss 
cone  corridor  is  about  850  km  long  (based  on  the  geometry  of  Fig.  1. 

See  also  Zlnn,  et.  al.) , and  that  500  km  is  needed  for  the  development 
of  the  turbulent  whistlers.  The  number  of  0+  ions  picked  up  by  the 
debris  ions  is  estimated  to  be  the  total  number  of  air  ions  in  the 
remaining  350  km  of  the  loss  cone  corridor.  Therefore,  we  have; 


Nf  (loss  cone) 


(0.35  x 108  cm) 
0.35  x 1013  A 


(A  cm^) 


(105) 


where  A is  the  cross  sectional  area  of  the  loss  cone  exit  and  where  we 
have  taken  nQ+  ^ 10-*  cm"^.  Assuming  that  A is  1%  of  the  surface  area 
of  the  parent  debris  bubble,  i.e.  A~.01  (4tiR2)~  1.13  x lO1^  cm2 

(for  R ~ 300  km),  we  have  that  Nf  (loss  cone)  ^ 0,4  x 1022  0+  ions. 
For  an  iron  bomb  weighing  106  gm,  Nfe~  1028,  only  a fraction  of  these 
particles  would  be  blown  into  the  loss  cone  exit.  It  is,  therefore, 
likely  that  air  ions  comprise  certainly  no  less  than  10%  of  the  plasma 
deposited  into  the  Van  Allen  belts  by  direct  injection  of  nuclear  debris 
from  a second  burst  into  the  loss  cone  exit  of  a pai  :nt  debris  bubble. 

It  is  possible,  however,  that  this  figure  could  be  much  larger,  perhaps 
100%  depending  on  the  fraction  of  the  second  bomb  debris  being  directly 
injected  into  the  loss  cone  exit,  the  length  of  the  loss  cone  corridor, 
and  the  actual  initial  debris  velocity. 
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The  length  of  the  loss  cone  corridor  in  which  shock-produced  accelera- 


Ik 


tion  of  ambient  ions  occurs  is  determined  ultimately  by  the  distance  the 
shock  travels  before  attenuating  to  a low  strength.  We  now  modify  our 
previous  estimate  by  including  this  attenuation.  In  this  calculation  we 
assume  the  debris  to  expand  outward  from  the  second  burst  into  the  existing 
loss  cone  corridors  of  the  parent  bubble.  We  know  that  the  debris -ambient 
ion  interaction  provides  an  environment  for  a two  stream  instability  that 
can  be  stationary  in  the  leading  edge  of  a propagating  shock  wave.  Con- 
sequently, we  assume  the  result  of  this  debris  motion  is  to  form  a shock  in 
the  loss  cone  corridor.  The  conservation  equations  for  this  collisionless 
shock  are  unaffected  by  the  magnetic  field  ( J_  to  shock  plane) , So  as  an 
estimate  of  the  shock  dynamics  we  can  consider  only  a gas  dynamic  shock  (if 
we  neglect  switch-on  type  shocks).  Thus  we  assume  that  a piston  of  mass/unit 
area  m with  a velocity  V0  moves  into  stationary  gas  at  conditions  p0,  pQ, 
uQ  ■ 0,  C0  “ CA  (CQ  ■ sonic  velocity).  At  t » 0,  the  piston  first  impacts 
the  ambient  gas  (Fig.  3). 


x “ 0 

I 


pressure  p(t)  front  with  pressure  pQ 


Figure  3 
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The  result  of  this  impact  Is  to  drive  a shock  wave  into  the  aabient 
gas  that  moves  with  speed  V8(t).  This  has  the  effect  of  decelerating 
the  piston  by  means  of  the  pressure  p(t)  exerted  on  the  piston  faces.  So 


mx  (t)  - -p(t) 
P 


(1) 


Now  we  assume  that  this  pressure  can  be  given  by 

2 


p(t)  2yM 


Y-l 


(2) 


P Y+1  Y+1 

o 

(pressure  ratio  across  a normal  shock) 

If  this  approximation  is  to  hold,  i.e.  that  p(t)  on  the  piston  face  is  the 
same  as  p behind  the  advancing  shock,  it  is  necessary  to  assume  that  the 
sound  speed  behind  the  shock  is  large,  i.e.  the  shock  is  very  strong. 
Therefore 


M » 1 

s 


(3) 


and 


p(t)  2 yM 


ar — 

Po  Y+l 


(4) 


Now  we  know  that  the  velocity  jump  across  a shock  is  given  by 


M 2 , i 2 

— = — • -U  — m 

.4.1'  1 w • 1 


Y+-1 


M. 


Y+1 


(5) 


where  [w]  = w2“  wi  is  the  velocity  Jump  (see  figure  4 ),  c^  is  the  sound 

speed  in  region  l and  is  the  upstream  (i.e.  unshockedl  Mach  number  in  the 

shock  front  rest  frame. 
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For  our  case  it  is  seen  that 

Cl  = C , w,  - Vs  , w2  = Vs-Vp>  Mx  = Ms(t)  and  [v]  = Vp(t) 
and  so  from  (5)  we  find 


V (t) 

_E 

c 

o 


Using  eq.  (7)  in  eq.  (2)  we  obtain 


Y+1  2 

P(t:)  * P0(  — )Vp  (t) 


where  we  used 


From  eqs . (1),  (8),  and  the  fact  that 


dV£> 

dt 


(10) 


xp(t)  = 


we  arrive  at  a nonlinear  differential  equation  for  the  piston  velocity 

yo 


(ID 


Equation  (10)  can  be  easily  solved  by  using  separability  and  the  initial  con 

dition  V (t*0)  = V . The  solution  of  eq.  (11)  is 

p o 


yt) 


(12) 


Now  we  assume  that  the  shock  will  stop  picking  up  particles  when  the 

shock  mach  number  M (t)  reaches  a mach  number  M-  ~ 2.76.  This  will  occur 
s 1 

at  a time  tf.  From  eqs.  (7)  and  (12)  we  find  that 


fcf  = 


m 


m 


(13) 


PoCoMf  P0v0  (Y+l)  P0CoMf 


where  eq.  (3)  was  used. 

From  this  attenuation  time  tf  one  can  obtain  the  effective  loss-cone 
corridor  length  over  which  pickup  occurs.  This  is  given  by 

r‘* 

fcf  » \ Vg(t)dt  (14) 

'o 


From  eqs.  (7)  and  (12)  we  obtain 
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m 

- — in 
f P 


V 1 y + 1 

1 + — — 

ca"j  2 


Noting  that  i fpQ  is  the  mass  pickup  per  unit  area  mf  and  that  &n 
is  of  order  unity,  we  find 


(15) 


i+ZalLli 

CAMf  2 


mf  5 m * mass  per  unit  area  of  the  piston  (16) 


Thus  the  debris  piston  is  predicted  to  pick  up  its  mass  in  air  ions. 

Before  examining  some  typical  numbers,  we  note  that  this  model  describes 
an  attenuation  mechanism  which  is  quite  different  from  a momentum  conserva- 
tion or  "snowplow"  type  model.  For  purposes  of  comparison  we  shall  present 
the  snowplow  idea  and  then  calculate  some  relevant  quantities  using  both 
models . 

For  a momentum  conserving  system  of  piston  plus  picked-up  air,  we  can 

write 


Vp(c) 


mv 


m+m  (t) 
a 


where  m (t)  is  the  mass  of  air  pickup  per  unit  area  at  time  t.  If  we  use  a 


strong  shock  approximation,  i.e.  eq.  (7),  we  obtain 

Y + 1 ® Vo 

m9w  - — „ 7 „s<t>  Tk 

When  pickup  ceases  Mg(t)  = Mf,  mfl(t)  = mf  and  we  find 


■ tn 


y + i 


- i 


2 MfC0 


y + i v 


m 


2 MfCA 


If  we  now  compare  the  two  models  we  find 


(18) 


(19) 
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m^  (Gas  Dynamic  Slowdown) 


(Snowplow) 


From  (20)  it  is  clear  that  m. 


(Snowplow) 


gives  the  higher  pickup  (we  shall 


see  it  is  typically  four  times  larger  for  a starfish- like  environment). 


>ical  Numerical  Results 


We  make  the  assumption  that  the  loss  cone  exit  comprises  1%  of  the 

6 2 

parent  bubble  area  A^  (—1.13  x 10  km*-)  and  so  total  debris  piston  mass  M 


(.01)^  = (.01)  (10J  kg)  =*  10  kg 


Assuming 


V = 1000  km/sec 
o 


(23a) 


C.  = 160  km/sec 
Ao 


(23b) 


Y * 5/3 


(23c) 


we  find  that  the  total  number  N of  oxygen  ions,  0 , picked  up  will  be 

m 

Hf  <G*S  Dy“"lc  S1°"d‘”"°  , .Ottj,  -^+  - .4  x 1027  particles  (24) 

^(Snowplow)  . _6  x l027  ^ (25) 

If  Vq  were  2000  km/sec,  estimate  (24)  would  not  change  significantly  due  to 
the  0(1)  property  of  the  in  function  (see  eq.  (15))while 


(Snowplow) 


= 1.65  x 10  particles 
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IV.  COMPARISON  OF  BURST  FLUXES  WITH  MAGNETOSPHERIC  SUBSTORM  FLUXES 


[ 

In  the  past  few  years  there  has  been  considerable  concern  over  the 
phenomena  of  synchronous  orbit  satellites  charging  to  high  potentials 
as  a result  of  magnetospheric  substorms  (Rosen,  1975).  These  substorms 
consist  of  the  injection  of  high  energy  plasma,  with  characteristics 
as  shown  in  Figure  5,  from  the  earth's  magnetotail  into  the  region  of 
synchronous  orbit.  This  injection  process  is  shown  in  Figure  6.  Those 
portions  of  a satellite  subject  to  the  high  energy  plasma  will  charge  to 
a potential  several  times  the  electron  energy  while  other  portions  of  the 
satellite  will  remain  at  ground  potential.  Potentials  near  ground  are 
maintained  by  photoelectron  emission  from  illuminated  surfaces  on  the 
spacecraft  or  by  contact  with  the  ambient  low  energy  plasma. 

During  eclipse  photoelectron  emission  disappears  and  during  a sub- 
storm the  ambient  low  energy  plasma  flux  is  strongly  dominated  by  the 
injected  high  energy  plasma.  The  most  damaging  discharges  as  a result 
of  spacecraft  charging  probably  occur  between  shadowed  spacecraft  com- 
ponents influenced  by  substorm  plasma  and  illuminated  components  at  ground 
potential.  When  the  discharge  passes  through  electrical  circuitry  between 
the  components,  damage  can  result.  Electromagnetic  interference  can  also 
result  from  surface  discharges  and  considerable  surface  deterioration 
can  be  caused  by  arcing. 

The  following  discussion  will  be  an  assessment  of  the  possible  space- 
craft charging  eiTocts  which  can  result  from  the  large  scale  transport 
of  ionospheric  plasma  to  synchronous  orbit  by  nuclear  bursts.  The  plasma 
instability  explored  in  the  previous  sections  of  this  report  demonstrates 
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Figure  5 Typical  Plasma  Electron  Characteristics 
as  Measured  in  the  Plasma  Sheet.  These  are  the 
electron  properties  attributed  to  the  midnight  to  dawn 
portion  of  synchronous  orbit  during  a substorm 
[After  Vasyliunas  (1968)] 


Figure  t Average  Particle  Motions  for  Various 
Energies  and  Particle  Types  during  a Substorm 
Plasma  Injection  from  the  Magnetotail  (After 
DeForest  and  Mcllwain  (1971)] 
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a mechanism  for  structuring  a collisionless  shock  wave.  This  mechanism 
will  operate  during  a high  altitude  nuclear  burst  as  the  expanding 
nuclear  debris  passes  over  the  stationary  ionospheric  plasma.  Through 
the  interaction  of  large  amplitude  whistler  mode  waves,  plasma  will  be 
picked  up  by  the  collisionless  shock  and  accelerated  into  the  magneto- 
sphere. 

As  the  large  mass  of  debris  and  plasma  propagates  along  the  earth's 
magnetic  field  lines,  the  plasma  mass  will  tend  to  expand  as  the  field 
lines  expand  at  high  altitudes.  This  is  illustrated  in  Figure  7. 

Using  nuclear  burst  parameters  for  a Starfish-like  (nominally  1.5 
megaton)  burst  (Zinn,  1966)  the  previous  calculations  have  conservatively 
estimated  that  near  1x10^  0+  ions  would  be  carried  by  the  shock  when 
17.  of  the  total  debris  exits  the  loss  cone.  Depending  on  orientation, 
mass,  and  shape  of  the  bomb  casing,  this  number  could  be  considerably 
higher. 


Magnetic  Field  Lines 


Synchronous  Orbit  Satellite 


Figure  7.  Nuclear  Debris  forms  collisionless  shork  wnich  transports 
Ionospheric  Plasma  to  higher  altitudes. 
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The  original  objective  of  this  research  was  to  study  the  injection 
of  plasma  into  the  magnetosphere  from  multiple  nuclear  bursts  and  their 
interaction.  The  greatest  yield  of  field  aligned  debris  fluxes  will 
probably  come  from  the  detonation  of  a nuclear  burst  inside  the  already 
formed  debris  bubble  from  a previous  burst.  If  the  bomb  casing  is 
cylindrical  and  with  its  axis  oriented  along  the  magnetic  field  lines, 
near  507,  of  the  debris  could  exit  the  loss  cone.  At  higher  altitudes  and 
at  night,  there  would  be  a great  proportion  of  the  ionospheric  plasma 

in  the  form  of  H+  with  a mass  of  one-sixteenth  that  of  0+.  The  debris 
shock  would  be  able  to  accumulate  significantly  more  of  this  plasma 
prior  to  turnoff  of  the  instability.  A conservative  worst  case  estimate 
of  the  total  ion  pickup  might  be  lxl 0^  ions,  with  more  than  an  order 
of  magnitude  contributed  by  the  proper  multi-burst  scenario  and  less  than 
an  order  of  magnitude  contributed  by  H+  pickup.  Plasma  neutrality  and 
electron  mobility  will  quickly  insure  that  the  p lasma  contains  an  equal 
number  of  electrons. 

The  total  propagation  distance  for  the  plasma  mass  from  the  burst 
point  just  above  the  earth's  surface  to  synchronous  orbit  is  approximately 
5x10^  km.  If  the  shock  is  initially  propagating  at  2000  km/sec  and  the 
Instability  shuts  off  with  a speed  of  1000  km/sec  and  free  streams,  it 
will  take  greater  than  50  sec  to  reach  synchronous  orbit.  This  time  is 
many  ion  cyclot.on  periods  at  all  altitudes  and  sufficient  to  have 
thermalization  of  the  complete  plasma  mass  at  near  the  ion  temperature 
as  the  mass  slows  and  diffuses  (Biskamp,  1973). 
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The  final  shock  velocity  after  instability  cutoff  and  free  streaming 
(an  estimated  500  km/sec  for  higher  altitudes)  corresponds  to  an  0+ 
energy  of  approximately  20  keV.  Little  of  this  energy  would  be  lost  in 
thermalizing  the  mobile  electrons  to  an  equivalent  energy. 

The  shock  structure  most  likely  forma  with  a spatial  extent  as  found 
by  (Zlnn,  1966)  and  as  mentioned  previously  of  500  km  along  the  field 
lines  and  of  30  km  radius  perpendicular  to  the  field  lines  for  the  loss 
cone.  Over  North  America  the  field  lines  for  synchronous  orbit  intersect 
the  earth  near  a magnetic  latitude  of  67°N  and  the  loss  cone  spreads  over 
one  degree  of  latitude. 

One  degree  of  latitude  at  67°N  will  span  an  L shell  difference 
of  .5  (L  ■ 6.3  to  L ■ 6.8).  At  synchronous  orbit  this  difference 
corresponds  to  a circular  area  of  expansion  for  the  loss  cone  of 
A^  ■ 8x10^  km2.  While  the  light  masses  of  the  electrons  will  allow  rapid 

thermallzatlon,  the  heavy  masses  of  oxygen  and  particular  the  iron 

debris,  will  constrain  their  initial  motion  closely  to  the  field  lines 
and  slow  crossline  diffusion.  The  outward  diffusion  of  plasma  caused 
by  plasma  pressure  and  the  diverging  field  lines  will  increase  the  radial 

scale  perpendicular  to  the  field  line  by  a factor  of  about  30,  and 

diffusion  along  the  field  lines  should  be  of  an  equivalent  degree. 

This  evaluation  of  the  spacecraft  charging  environment  from 

multiple  nuclear  bursts  la  therefore  estimated  to  be  a total  of  10 

electrons  of  20  keV  energy.  These  electrons  would  be  spread  over  an 

6 2 

area  of  « 8 x 10  km  at  synchronous  orbit.  They  would  be  in  a 

4 

shock  piston  of  length  L - 30  x 500  km  ■ 1.5  x 10  km  that  would 

have  a duration  at  a synchronous  satellite  of  t • 1.5  x 10^  km/500  km/sec  - 

30  sec. 
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Based  on  Che  previous  assumptions  as  Co  Che  nature  of  plasma  dep- 
osition by  multiple  nuclear  bursts  It  is  possible  to  compare  these 
fluxes  with  observed  magneto spheric  substorm  electron  fluxes.  A total 

of  10^  electrons  spread  through  a volume  of  V ■ A^  x L - 1.2  x lO1^3 

2 3 

gives  a density  of  approximately  8 x 10  electrons/cm  . This  density 

12 

at  an  energy  of  20  keV  translates  to  an  omnidirectional  flux  of  6.7  x 10 
2 

electrono/cm  -sec.  These  fluxes  are  significantly  greater  than  those 
encountered  during  a magnetospheric  substorm,  as  shown  in  Figure  8. 

While  the  duration  of  the  nuclear  injection  is  much  shorter  than  a 
typical  subatorm,  studies  (Rothwell,  et  al,  1976)  indicate  the  charging 
process  Cakes  only  fractions  of  a second. 

Figure  8 The  Omnidirectional 
Flux  of  Electrons  with  Energies 
(650  eV  < E < 53  keV)  during  a 
Substorm,  with  a Standard  Reference 
Day  for  Comparison.  [After  Sharp 
and  Johnson  (197  2)] 

In  both  the  single  and  multiple  nuclear  burst  environments  there 
appears  to  exist  the  potential  for  spacecraft  charging  effects  at  syn- 
chronous orbits.  The  calculated  fluxes,  densities,  and  energies  of 
the  Injected  electrons  are  all  greater  than  those  encountered  during 
substorms  and  which  have  been  observed  to  cause  spacecraft  charging. 

These  calculations  contain  many  approximations,  but  until  further 
research  can  refine  the  results  it  may  be  expected  that  synchronous 
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orbit  satellites  under  certain  nuclear  burst  conditions  will  find 


themselves  subject  to  a short,  but  Intense,  period  of  spacecraft 
charging.  Potentials  In  the  tens  of  kilovolts  for  periods  of  tens 
of  seconds  are  suggested.  The  resulting  transient  charging  and  arc 
discharging  should  surely  be  a part  of  the  design  criteria  of  any 
spacecraft  expected  to  survive  a situation  where  high  altitude 
nuclear  bursts  are  Involved. 
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. DISCUSSION  AND  RECOMMENDATIONS 


The  results  of  this  research  have  shown  that  a multiple  nuclear  burst 
can  provide  the  necessary  environment  to  excite  an  ion  cyclotron  beam  mode- 
whistler  mode  two  stream  instability  that  is  stationary  in  the  rest  frame 
of  weak  parallel  shock  waves.  The  expansion  of  debris  from  such  bursts 
can  therefore  produce  a collisionless  shock  wave  in  the  loss  cone  corridors 
of  an  existing  nuclear  debris  bubble.  Estimates  of  the  resulting  electron 
fluxes  along  the  distended  geomagnetic  field  lines  into  the  magnetosphere 
indicate  that  this  mechanism  could  produce  20  kev  omnidirectional  electron 
fluxes  of  approximately  6.7  x 1012  cm“2  s_1  for  durations  as  high  as  30  s. 
These  fluxes  are  two  orders  of  magnitude  larger  than  those  due  to  magneto- 
spheric  substorms  and  occur  for  times  one  order  of  magnitude  longer  than 
that  required  for  charging  of  synchronous  orbit  spacecraft.  It  would  thus 
appear  that  multiple  nuclear  bursts  could  be  a significant  source  of  ener- 
getic particle  fluxes  into  the  magnetosphere  and  the  cause  of  consequent 
charging  of  synchronous  orbit  spacecraft. 

Since  these  conclusions  are  based  on  a linear  dispersion  theory  of  the 
relevant  two  stream  instability  and  relatively  simple  calculations  of  par- 
ticle fluxes  caused  by  the  passage  of  the  resulting  collisionless  shock 
wave,  two  lines  of  further  study  are  suggested: 

(i)  Continue  the  study  of  the  hot  ion  cyclotron  beam-whistler  into 
the  quasi-linear  stage  of  its  development  with  a calculation  of 
the  saturation  level  of  the  turbulence. 

(ii)  Improve  calculations  of  particle  fluxes  by  incorporating  the 
results  of  the  nonlinear  study  and  by  modelling  the  particle 
pick-up  mechanism  more  realistically. 
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Crossfield  Magnetosonic  Two  Stream  Instability1 
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In  (his  article,  a detailed  analysis  was  made  on  the  plasma  magnetosonic  instabilities  whose 
propagation  is  perpendicular  to  the  magnetic  held  and  whose  two  streams  cross  the  magnetic  field 
during  the  time  of  development  of  the  instability 

We  have  also  studied  the  role  of  the  magnetosonic  instabilities  in  the  structure  of  perpendicular 
shock  waves.  Since  these  instabilities  are  not  stationary  in  the  shock  interior,  one  can  ultimately 
deduce  that  they  do  not  play  a pnncipal  role  in  the  structure. 


Dans  cet  article,  une  analyse  detaillee  aete  faite  sur  les  instabil.tes  magnetosonores  de  plasma 
dont  la  propagation  est  perpendiculaire  au  champs  magnetique  et  donl  les  deux  ecoulements 
traversent  le  champs  magnetique  pendant  le  temps  de  developpement  de  I'instahilite. 

Nous  avons  aussi  etudie  la  role  des  instabilities  magnetosonores  dans  la  structure  des  ondes  de 
choc  perpendiculaires.  Comme  ces  instabilities  ne  sont  pas  stationnaires  a l interieur  du  choc,  on 
peut en deduire  finalment  qu  elles  nejouent  pas  une  role  pnnctpale  dans  la  structure. 

Can.  J.  Phys.,53.  1022(1975) 

I.  Introduction 

The  counterstreaming  of  two  ion  beams  across 
a perpendicular  magnetic  field  in  the  presence  of 
a stationary  electron  background  gives  rise  to  a 
spectrum  of  unstable  growth  modes  that  were 
thought  to  provide  the  necessary  turbulent  dissi- 
pation to  structure  weak  collisionless  shock 
waves  (Papadopoulos  et  at.  1971).  In  this  note 
a complete  study  of  this  magnetosonic  instability 
is  presented,  differing  from  the  linear  analysis  of 
Papadopoulos  et  al.  (1971)  in  that  (i)  the  com- 
plete magnetosonic  dispersion  curve  is  used 
(resulting  in  a wider  range  of  relevant  drift  Mach 
numbers),  (//)  the  weak  beam  - magnetosonic 
mode  interaction  is  considered  as  well  as  the 
equidense  beam  theory,  and  (Hi)  pertinent  cal- 

[1]  « 1 « k*v*iaji  nj  « m2  « Qj 

where  is  the  ion  (electron)  cyclotron  frequency.  Furthermore,  the  parameter  ordering  is 

12}  a.  « to  « ka,  « « copf;  to  « Ac;  Qct  « |a>  - AK,|  (s  = 1,  2) 

where  wpl(t)  is  the  ion  (electron)  plasma  frequency  and  at  is  the  electron  thermal  velocity.  The 
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culations  have  been  performed  for  the  t*  folding 
growth  lengths  of  unstable  modes  br  th  at  the 
leading  edge  and  density  midpoint  of  weak 
collisionless  shock  waves. 


II.  Linear  Analysis 

In  this  analysis,  the  two  ion  beams,  with  den- 
sities rt,  and  velocities  Vx  (s  = 1,  2),  are  assumed 
to  be  unaffected  by  the  magnetic  field  during  the 
full  development  of  the  instability  while  the 
electrons  are  s rongly  bound  to  the  field  lines, 
so  that  the  unperturbed  system  consists  of  two 
uniform  beams  with  constant  drift.  The  wave 
number  and  complex  oscillation  frequency  there- 
fore must  satisfy 


I 


.. 


1 
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perturbed  plasma  conductivity  a is  found  from  the  linearized,  Fourier  transformed  fluid  equations 
for  ions  and  electrons  in  conjunction  with  Faraday’s  law.  Used  in  the  dispersion  relation, 

e„  = ex,2(D2//c2c2  (£  = 1 + 4jti'S/co) 

where  x and  y directions  are  taken  along  F,||A  and  B0  respectively,  this  in  turn  leads  to  the  dis- 
persion relation 


In  the  electron  rest  frame,  the  charge  and  current  neutrality  condition  of  the  unperturbed  system 
can  be  incorporated  into  the  single  equation,  (1  - q)K,  + qF2  = 0 where  q = n2/n,  is  the  beam 
concentration.  Then  upon  introducing  into  [3]  the  relative  drift  velocity  U = V1  — Vt  and  the 
effective  ‘plasma  frequency’, 

«„(*)  = [ftcAi/O  + o>,k 2/*2c2)],/2 

one  obtains, 

[4]  (i  - n)o>02(*Xa>  + n kuy1  + n<u02(*)[<u  - (i  - r\)kuy2  = 1 

In  this  convenient  form  both  the  weak  (q  « 1)  and  equidense  (q  = 1/2)  beam  theories  may  be 
analyzed. 

Consider  first  q « 1.  Setting  q = 0 in  [4]  leads  to  the  usual  magnetosonic  dispersion  relation, 

to  = a>c(k)  = kCA(  1 + k2c2/ CDp,2)" 1/2 

where  CA  is  the  local  Alfven  speed.  Possible  strong  interactions  may  be  expected  in  the  vicinity 
of  the  intersection  of  io0(k)  with  the  beam  dispersion  curve  tob  = kU  — kV1.  Since  the  slope  of 
<o(k)  satisfies  the  inequality  to '{k)  < o)'(0)  = CA,  it  is  clear  that  such  interactions  can  occur  only 
if  = U/CA  ^ 1.  A perturbation  analysis  of  the  full  dispersion  relation  for  q « 1 in  the  vicinity 
of  the  intersection  shows  that  the  spectrum  of  unstable  modes  is  given  by 

0 < k1  s?  {My1  - I K.V  = (Mt1  « I) 

with  maximum  growth  rate  occurring  at  k = kc  and  given  by 

[5]  y = Im  (to)  = q‘^3I/22'4/J[(l  - /Wd2)ncin„l‘/J 

The  equidense  beam  case  can  be  analyzed  by  setting  q = 1/2  and  solving  the  resulting  quartic 
equation  for  w2(£).  It  is  found  that  for  the  k spectrum  given  by 

0 < k1  < (4A/d  - 2 - I )opc2/e2  {Ht  = U/cA) 

w is  pure  imaginary  with  growth  rate 

[6]  y{k)  = + 2Md2(l  + kVl^V11  - 1 - (1/2)M„2(1  + *2c2/g>p,2)}  1/2 


This  is  equivalent  to  the  corresponding  result  of 
Papadopouios  et  at.  (1971)  only  in  the  limit 
k2c2  » o)pt2.  The  mode  with  maximum  growth 
is  then  found  from  dw/d/c  = 0. 

III.  Shock  Structure  and  Discussion 

To  examine  the  implications  of  this  magneto- 
sonic  instability  for  collisionless  shock  structure, 
a shock  layer  is  modeled  as  a region  of  inter- 
penetrating upstream  (unshocked)  and  down- 


stream (shocked)  ions  (e.g.,  Golden  el  al.  1973). 
In  the  reference  frame  of  the  leading  edge  of  the 
shock  the  ion  velocity  distribution  takes  on  the 
bimodal  Mott-Smith  form 

/,(*,  V)  = n„(x)8(F  - yj  + r.d(x)S(F  - Vt) 
and  for  the  electrons 

/,(*,  V)  = "«(x)6(  V - yt(x)) 

In  these  distributions  Pu  and  Vt  are  the  average 
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upstream  and  downstream  ion  velocities  and 
/!„(*)  and  «d(x)  are  the  respective  ion  densities 
through  the  shock.  The  corresponding  electron 
quantities  are  determined  by  the  charge  and 
current  neutrality  conditions  as 

n.(x)  = nu(x)  + nd(  x) 

K(x)  = [1  - t)(*)]Ku  + q(x)Fd 

where  nW  = Hi(x)/nt(x). 

The  results  of  the  preceding  section  may  now 
be  applied  to  the  shock  structure  problem  by 
Doppler  shifting  the  frequency  w according  to 
to  = to,  — fcF.and  using  the  Rankine-Hugoniot 
relations  for  a magnetic  shockwave.  The  growth 
rate  of  the  weak  beam  theory  may  be  used  to 
examine  the  behavior  of  unstable  modes  in  the 
leading  edge  of  the  shock.  If  this  instability  is 
to  play  a role  in  structuring  a weak  collisionless 
shock,  then  unstable  modes  should  either  have 
zero  group  velocity  in  the  frame  of  the  leading 
edge  or  should  propagate  slowly  downstream. 
Although  the  modes  are  found  to  propagate  into 
shock,  even  the  maximum  growth  rates  are  suf- 
ficiently small  that  the  modes  must  travel  an 
appreciable  fraction  of  the  shock  thickness  before 
experiencing  an  e fold  increase  in  amplitude.  A 
typical  calculation  of  the  e folding  distance 
A = y~idu>Jdk  for  the  maximum  growth  mode 
near  the  leading  edge  of  an  air  plasma  shock 
shows  that  A = 7.3 c/u>peo  (subscript  'o’  refers  to 
leading  edge  conditions)  for  an  upstream  Alfven 
Mach  number  MK  = 2l/I  and  yt  = cp/cv  = 2.  In 
this  case  a choice  of  q = (m./w*)1'2  is  simul- 
taneously small  enough  to  insure  validity  of  the 
perturbation  expansion  yet  not  too  small  to 
violate  the  last  of  inequalities  [2],  Since  the  shock 
thickness  L,  ~ 12c/a>peo  (Golden  et  al.  1961),  it 
is  unlikely  that  this  leading  edge  mode  contri- 
butes to  the  structure  of  the  shock. 

It  is  also  useful  to  consider  the  e folding  growth 
distances  for  unstable  modes  generated  in  and 
propagating  into  the  midshock  region  (q  = 1/2). 
In  the  frame  of  the  shock  leading  edge  it  is  found 
from  [4]  that  d<t>Jdk  = Fu  so  that  A = ■/'’(£)  Fu. 


In  particular  for  an  air  plasma  shock  of 
Ma  = 21/J  and  yt  = 2 two  typical  conclusions 
may  be  reached:  (i)  for  the  leading  edge  mode  of 
maximum  growth,  assumed  to  reach  the  shock 
midpoint  while  still  in  the  linear  stage,  the  e 
folding  length  is  reduced  only  to  4.4c/(0pto  so  that 
the  variation  of  A through  the  shock  is  not  suffi- 
ciently strong  to  achieve  much  shorter  e folding 
lengths,  and  (it)  the  maximum  growth  mode  at 
q = 1/2  produces  a minimum  A of  3.6(c/<opttl)  so 
that  even  this  mode  can  only  grow  1^-2  e folds 
as  it  propagates  from  the  midshock  region  to  the 
trailing  edge.  It  should  be  noted  that  these  results 
are  typical  and  relatively  insensitive  to  the  choice 
of  yt  and  to  the  shock  strength. 

IV.  Conclusions 

The  results  of  this  study  have  been  a complete 
analysis  of  the  magnetosonic  two  stream  insta- 
bility including  growth  rates,  limits  of  validity, 
and  unstable  spectra  for  both  the  weak  beam- 
magnetosonic  mode  interaction  and  the  equi- 
dense  countrrstreaming  beam  instability.  Fur- 
ther, it  has  been  found  that  this  instability  is 
probably  not  an  important  source  of  turbulent 
dissipation  fer  the  structure  of  weak  shocks.  On 
the  other  hand,  for  shocks  where  the  initial  steep- 
ening is  broadened  by  classical  dissipation,  the 
instability  might  play  a role  in  broadening  the 
trailing  edge.  As  a consequence  it  is  now  believed 
that  perpendicular,  collisionless  shocks  are 
structured  by  turbulence  originating  from  the 
so  called  current  driven  modified  two  stream 
instability  (Ott  et  al.  1972). 
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APPENDIX  B 


The  attached  Physics  Letter  reprint  summarizes  some  initial 
aspects  of  our  research  on  the  application  of  streaming  plasma 
instabilities  to  parallel  shock  wave  structures.  A complete 
discussion  of  this  work  is  given  in  Sec.  B of  this  report. 
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A parallel  shock  structure  is  modeled  as  regions  of  interpenetrating  streams  of  cold  unshocked  ions  and  hot  shocked 
ions.  Our  linear  dispersion  theory  predicts  that  unstable  whistlers  can  stand  at  the  leading  edge  of  weak  and  interme- 
diate strength  shocks. 


Recently  [1],  it  was  proposed  that  unstable  beam-whistler  wave  interactions  formed  the  turbulent  structure  of 
shock  waves  propagating  along  the  magnetic  field  (parallel  shocks)  in  collisionless  plasmas.  It  was  assumed  that 
the  electrons  could  be  modeled  as  warm  fluid  and  the  ions  as  a Mott  Smith  superposition  of  the  unshocked  and 
shocked  ion  flows:  for  mathematical  simplicity,  these  ion  flows  were  modeled  as  cold  monoergetic  beams.  Ana- 
lysis of  the  ensuing  linear  dispersion  relation  then  revealed  the  existence  of  unstable  modes  standing  near  the 
leading  edge  of  the  shock  for  Alfven  Mach  numbers  > 5.5.  This  suggests  that  such  modes  are  most  likely  to 
grow  to  sufficiently  large  amplitude  to  scatter  incoming  ions  and  create  the  required  dissipation  for  intermediate 
strength  shocks.  It  was  recognized,  however,  that  the  assumption  of  cold  ion  streams  was  a serious  limitation  in 
the  theory  since  (i)  the  shocked  downstream  tons  are,  in  any  case,  thermalized,  and  (ii)  the  upstream  ions  must 
be  sufficiently  warm  (rs  = \Zyn(Te  + 7'j)/mi  > CA  = B0l\/4trpu ; B0  ~ constant  magnetic  field,  pu  = upstream  mass 
density)  to  preclude  the  possibility  of  “switch-on"  shocks  (rotation  of  the  magnetic  field  across  the  shock  layer). 

In  this  letter,  we  pa  tially  relax  the  cold  ion  stream  requirement  of  the  previous  theory  by  more  realistically 
allowing  the  downstream  ions  to  be  Maxwellian  with  final  temperature  Td  and  mean  velocity  Fd  dictated  by  the 
Rankine-Hugoniot  relation.  Then  following  the  the  Mott-Smith  formalism,  the  ions,  in  the  rest  frame  of  the  shock 
front,  are  modeled  as  the  superposition  of  two  streams:  hot  downsteam  particles  moving  with  relative  velocity 
(Kd  - Vu ) through  the  cold  upstream  particles.  The  appropriate  velocity  distribution  function  for  the  ions  at  the 
station  z inside  the  shock  is  therefore  given  by 


/i(z,u)  = wu(z)5(u-  Fu)  + 


”d<z> 


exp 


- (u-Fd)21 

C 2 J’ 

<-d 


0) 


where  Cd  = \/2 The  warm  background  electrons  move  with  mean  velocity  Fe(z)  = ( 1 - tj(z)]  Fd  + rrtz)Ft, , 
(t)  = nd/(nu  + dictated  by  the  requirements  of  local  charge  and  current  neutrality.  Only  weak  beam-plasma 
interactions  (rj  < 1)  arc  considered  here  (applicable  therefore  to  the  leading  edge  of  the  shock  layer)  and  consider- 
ations of  arbitrary  hot  beam  strengths  are  deferred  to  a later  paper. 

For  the  case  of  parallel  propagation  (A||Z?0||z-axis),  the  linear  dispersion  relation  is,  in  the  electron  rest  frame. 


exx(kio)  ± itxy(ku)  = (kc/to)2  . 


(2) 


The  dielectric  components  in  (2)  are  calculated  by  addition  of  the  ion  and  electron  polarizabilities.  These  well-known 
polarizabilities  are  formulated  explicitly  in  terms  of  the  (unperturbed)  ion  and  electron  distribution  functions  and 
therefore  take  account  of  the  bimodal  nature  of  the  ions  (eq.  (1)  in  the  shock  reference  frame)  and  the  locally 
Maxwellian  behavior  of  the  electrons.  We  make  the  following  low  frequency  approximations:  (i)  the  vacuum  dis- 
placement current  is  negligibly  small  (w  < kc)\  (ii)  the  electrons  have  zero  mass  (Icol  < nce  = eB0lmec)  so  that 
they  can  be  modeled  as  warm  fluid  (krLe  < 1).  The  dispersion  relation  in  the  reference  frame  of  the  shock  front 
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then  ultimately  becomes 

d u exp[-Q</Cd)2]  + k2C\(z)  kVt(z) 

^-kV^tlJ-kU  J2ci  fig 

where  u>,  = cj  + k • Ve  is  the  frequency  in  the  shock  frame.  Note  that  the  cold  dispersion  relation  (5)  yf  ref.  [ 1 ] is 
recovered  from  (3)  by  letting  C'd  -*  0 and  transforming  back  to  the  electron  frame.  The  solution  of  eq.  (3),  cj,  = 
cjt(k),  should  provide  the  spectrum  of  unstable  growth  modes.  Only  those  growth  modes  which  have  zero  group 
velocity  (in  the  shock  front  rest  frame)  are  of  interest,  since  only  these  have  sufficient  time  to  grow  to  large  am- 
plitude to  scatter  tncoming  particles. 

In  the  weak  beam  (17  < 1)  approximation,  we  set  cjs=  v + i-y,  7 < |p|.  The  solution  of  (3)  is  then  readily  found 
to  be: 

K*;  yu)  = kvu  + ^2c^/(2fici)  - kCxs/iV(kcj2nci)2 , . 

7(*;^u)  Vffflci  (p-JtKu+flJ2  r (p-^^+flg)2-! 

flci  " l*ICd  ^ {v-kVJ(.V-kV^2^)^P  L k2C\  j' 

Strong  interaction  evidently  corresponds  to  the  intersection  of  the  v - kKd  + £2^  = 0 hot  ion  beam  line  with  the 

cold  whistler  dispersion  curve,  eq.  (4),  in  the  third  quadrant  of  the  (cus,  fc)-plane.  Here  the  resulting  7max  varies  j 

directly  with  the  strength  •}  of  the  hot  beam.  In  the  previous  theory  [1  J,  where  the  shocked  ions  were  modeled  1 

as  a cold  monoenergetic  beam,  the  corresponding  weak  ion  beam-whistier  interaction  was  found  to  be  stronger, 

i.e.  7max  a >/n-  This  is  not  surprising  since,  in  the  present  theory,  the  whistler  wave  interacts  strongly  only  with  a 

small  percentage  of  the  hot  beam  ions,  whereas  in  the  previous  theory,  all  of  the  particles  in  the  monoenergetic  j 

ion  beam  participate  in  the  interaction. 

The  (rj  < 1)  zero  group  velocity  condition  bv(k,  Vu)ldk  = 0,  readily  obtained  from  (4),  reveals  that  a given 
shock  velocity  > CA  corresponds  to  a given  value  of  k < 0.  From  (4)  and  (S)  we  see  that  k < 0 is  suffi- 
cient for  >(*)  > 0,  so  that  whistlers  standing  at  the  leading  edge  grow  there  (this  instability  does  not,  for 
the  most  part,  correspond  to  an  intersection  between  the  hot  ion  beam  line  and  the  whistler  dispersion 
curve)  This  is  evidently  true  for  all  MA  > 1 up  to  the  limit  of  validity  of  our  me  = 0 model  (given  by  MA 
<(2/3)V"ij/me,  see  ref.  [1]).  Thus  the  turbulent  whistler  mode  structure  for  intermediate  shocks  (4fA  > 5.5) 
suggested  by  the  previous  cold  beam  theory  is  suggested  as  well  for  weak  shocks  (K  AfA  < 5.5)  by  the  present 
hot  beam  theory. 
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Equation  (3)  should  be  modified  to  read: 


n - i<*>l  — dr-nr-  * "(z)  ~=—  f du  c-(u/cd)z]  j- 

“s  - kVu  + -fc,  ) d 


si  . - ku 
ci 


K cf(z)  a)  kV  (z) 

= 1 + 2 _ _ + s 

i,  °«i  ' 


Equation  (5)  should  be  modified  to  read: 


Y(k;V„) 


(v  - kV  + Si  .)  (k V.  - v) 


u ci 


|k|Cd  (V  - kVu)(v  - kvu  + 2nci) 


(v  - kvd  + Qciy 


Page  252,  last  paragraph,  on  lines  2 and  3 instead  of  "From  (4)  and  (5), 
we  see  that  k < 0 is  sufficient  for  y(k)  > 0,  so  that  whistlers  standing..." 
should  be  written  "From  (5),  we  see  that  v(k;Vu)/Vd  < k < 0 is  sufficient 
for  y(k)  > 0,  so  that  certain  whistlers  standing  " 

Page  252,  last  paragraph,  on  line  5 instead  of  M.  >_  1 should  be  written 
Ma>2.76. 

Page  252,  last  paragraph,  second  line  from  bottom,  instead  of 
1 <.  _<  5.5  should  be  written  2.76  <_  £ 5.5. 

These  modifications  do  not  alter  the  main  statement  of  our  paper. 
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